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Foreword 


The problem of providing protection against the radiations generated by 
nuclear weapons is of much concern to the Federal Government and to the 
people of the United States. A great deal of effort and money is being spent 
currently on attempts to analyze existing structures for shielding properties, 
and plans are being carried forward to design and build special shelters of 
many kinds. 

The National Bureau of Standards has maintained a continuing research 
program directed to the development of engineering data and methods basic 
to the solution of these radiation shielding problems. Most of the work of the 
past few years on this project has been concerned with gamma rays from radio- 
active fallout. This Monograph summarizes both the methods and the data 
which have been generated as a result of research on the penetration of fallout 
radiations. 

Financial support for this project has come primarily from the Office of 
Civil and Defense Mobilization and from the Defense Atomic Support Agency 
of the Department of Defense. 


А. V. Astin, Director. 
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Structure Shielding Against Fallout Radiation From Nuclear 
Weapons 


L. V. Spencer 


The theory of structure shielding from fallout gamma radiation is developed to the point 


of applications to elementary structure types. 


Examples discussed in the text include the 


density interface, foxhole, shielded foxhole or basement, light superstructure, vertical wall, 


blockhouse, vents, compartmentalization effects, and mazes. 


A large number of engineer- 


ing charts and graphs are presented for engineering calculations, including many obtained 


from angular distributions of the exposure dose. 
This information has been obtained almost completely 


and for Со and С5137 sources. 


by machine calculations utilizing basic cross section data. 


Results are given for a fission spectrum, 


A number of sources of experi- 


mental data are mentioned, but detailed comparisons with experiment are not included. 


I. Introduction 


1. Purpose 


This Monograph is to assist scientists and engi- 
neers in the solution of problems of protection 
from ionizing radiation, particularly radiation 
from fallout. Important activities in which such 
problems arise include (a) the design of.new pro- 
tective structures, (b) analysis of existing struc- 
tures to determine the natuie of their shielding 
capabilities, and (c) assessing the possibilities for 
improvising shelter. 


2. Comments about Radiation Shielding 
Problems and Publications 


The analysis of structures for their protection 
against ionizing radiation represents a new field 
of engineering which has many similarities to 
illuminating engineering. Estimating radiation 
levels at different locations in a complex structure 
is comparable to the problem of determining the 
illumination levels in a similar stiucture on a 
cloudy day, but with all partitions and walls 
having varying degrees of transparency or trans- 
lucence rather than being opaque. The com- 
plexity of the problems is such that, despite 
intensive research, significant gaps in our knowl- 
edge still remain. 

Reactor shielding problems are rather different 
from the problems of shielding against nuclear 
weapons. The former involve a localized source 
and shield whose various aspects are readily 
available for study and possible modification. 
Further, both source and shield differ in many 
ways from the sources and the barrier configura- 
tions encountered in weapons shielding problems. 
For these and other reasons the different reactor 
shielding handbooks are not an adequate basis for 
the study of weapons shielding problems, though 
Е: do contain a great deal of useful information 

1—5]. 
The joint AEC, DOD, and OCDM publication, 


“Тһе Effects of Nuclear Weapons," gives a general , 


description of all aspects of nuclear explosions, 


*Figures in brackets indicate the literature references at the end of this 
Monograph. 


together with many tables and graphs relevant to 
radiation shielding [6]. However, the material 
included in the book on the subject of protection 
against fallout constitutes an introduction to the 
problems rather than a detailed methodology for 
their solution. 

A number of reports have been published in the 
past three years on shielding against fallout 
radiation [7-11]. The problems and the different 
approaches to their solution contained in these 
reports tend to have close similarities to each 
other and to the material in this monograph. 
Nevertheless, there are differences worth noting. 
The approach taken here is less empirical, relying 
on basic data derived from first principles rather 
than on field test experiments. This has the 
advantage of providing more different types of 
data, as well as a basis for different kinds of 
estimates. Enough experimental confirmations of 
the basic data exist to give confidence in its overall 
accuracy [12]; but further research can be expected 
to complement and modify this information in 
various ways.! 


3. Choice and Presentation of Material 


The contents and organization of this Mono- 
graph have been influenced by the fact that the 
analysis of structure shielding may properly be 
considered a branch of either nuclear or civil 
engineering. It is, however, a subject with which 
few people are familiar; hence it has seemed worth 
while to include general information about fallout 
radiation and its characteristic properties. Thus, 
Part II describes fallout in an elementary way and 
Part III gives properties of gamma radiation. 

Part IV introduces the formal approach to 
structure shielding analysis adopted іп this 
manuscript together with the basic parameters of 


1 Recently, OCDM has published a series of manuals, prepared by L. N. 
FitzSimons of that agency, and based on the methods and data presented in 
this monograph. See, e.g., ‘‘Fallout Shelter Surveys: Guide for Architects 
and Engineers," Report No. N P-10-2 (May 1960), and "ОСОМ Engineering 
Manual,” (Dec. 1, 1960). 


description and measurement. Part V isa dis- 
cussion of graphical data for use in applications. 
Part VI gives an analysis of elementary structures. 
Finally, Part VII presents a methodology for 
generalizing the calculations to a wide group of 
barrier shapes. 

It was originally intended to incorporate into 
this Monograph analyses of a number of complex 


structures in terms of the data presented here. 
However, this has been omitted for a number of 
reasons. Instead, we include an appendix which 
gives two additional sets of graphical data, one 
for Со and the other for Cs", for use in analyzing 
experiments. The data contained in the body of 
the manuscript are for radiation from fission 
products. 


II. Description of Fallout 


4. Fallout vs Initial Effects 


Whenever a nuclear bomb is exploded, a tre- 
mendous pressure wave, large quantities of heat 
and light, and a vast amount of radiation and 
radioactive material are produced concurrently 
|6, Chapters II and III]. Where these three 
effects are all important, their interrelations com- 
plicate nearly all prctection problems. For ех- 
ample, completely protective structures must be 
blast resistant as well as radiation resistant, and 
they must be built with due consideration for the 
possibility of a fire storm. Furthermore, the 
nature of blast and heat phenomena affect the 
radiation problem, e.g., through vertical drafts 
which alter the spatial distribution of radioactive 
material. 

Exposure to the initial radiations ceases about 
a minute after the explosion, by which time the 
radioactive material from the explosion is carried 
out of range by the upward drafts. Large 
amounts of earth and debris are also drawn up- 
ward, and the resulting mushroom-shaped cloud 
may rise 80,000 feet or more. ‘This cloud contains 
particles ranging in size from submicroscopic 
specks to grains as large as coarse sand. Radio- 
active residue from Ше bomb and its surroundings 
adheres to these particles. The general term for 
this particulate matter, as it is brought to earth 
again, is “fallout.” 


5. Regional Variations in Fallout 


Some of the fallout particles rise higher than 
others; and some are more buoyant than others, 
thus falling to earth less rapidly [6, Sections 
10.1-10.24]. Close to the explosion there occurs 
the “throwout’’—rocks blown out by the blast; 
while at the other extreme are particles so small 
and raised so high above the earth that they are in 
the stratosphere above the weather and are little 
influenced by gravitational forces. It takes years 
for these tiniest particles to settle to earth again, 
and they may come down anywhere on the earth. 
We shall not discuss this long-duration ‘werld- 
wide" fallout, nor shall we be concerned with the 
"throwout" [6]. 

А large fraction of the fallout material is big 
enough to come to earth within hours or days, yet 
small enough to be easily blown about by winds. 
Correspondingly, the distribution of particle 


200 


300 ROENTGENS 


190 


180 


170 


160 


150 


140 


130 


120 


IIO 


100 


90 


во 


то 


DISTANCE FROM 6.2. (MILES) 


60 


50 


40 


30 





го 10 О о 20 


DISTANCE FROM 6.2. (MILES) 


FIGURE 5.1. Idealized total (accumulated) dose contours 
from fallout in the first 36 hours after the high yield explo- 
sion at Bikini Atoll on March 1, 1954. (Ref. 6.) 


sizes and the regional wind structure determine 
the spatial distribution of the deposit. Sustamed 
strong winds may carry fallout in lethal amounts 
for hundreds of miles, whereas light winds may 
leave most of 1t much closer to the explosion. 

Figure 5.1 gives idealized contours of total fall- 
out exposure from a large (multi-megaton) ex- 
plosion, indieating how the fallout was carried by 
winds away from the point of explosion. 
Attempts to calculate such gross variations have 
had some success, at least to the point of being 
useful in estimations for broad planning purposes 
|6, Sections 9.559.138]. 

It is clear from figure 5.1 that intense fallout 
can oceur in wide areas far beyond blast and heat 
effects. This makes it reasonable to consider the 
fallout as a problem separate from blast and heat 
complications. In this Monograph, problems of 
protection from fallout radiation are discussed 
without reference to the effects of blast and heat; 
it should be remembered, however, that intense 
fallout also occurs in the proximity of the 
explosion. 


6. Local Variations in Fallout 


Just as the gross features of fallout spatial 
distributions are affected by the general wind 
pattern over a large area, so the fine features 
depend on the local winds. For example, fallout 
may collect more on one side of a building than 
on another, or blow down one set of streets more 
than another, more or less in the manner that the 
drifting of snow is determined by local variations 
in the surface winds. We do not have explicit 
data on the effects of such variations, and a 
thorough analysis of the problem is needed before 
a reliable estimation of the magnitude of these 
"micro-meteorological" effects can be established 
[13]. In the absence of specific information we 
assume for the present that fallout will be uniformly 
distributed horizontally over exposed surfaces. | 

А question of рагиещаг importance relates to 
the amount of fallout material which may remain 
suspended above the earth in trees and shrubbery. 
This has not been accurately determined. There 
is evidence that not more than a few percent of 
the material descending, e.g., upon a tree will 
remain on the limbs and leaves of the tree if the 
surfaces are dry [14-15]. But the possibility of 
weather conditions which will cause fallout to 
adhere even to vertical wall surfaces, and the 
sensitivity of basement protection to the amount 
of suspended source material, make this a major 
problem for future research. 


7. Types of Fallout Radiation 


Two types of radiation are emitted by radio- 
active fallout materials,? namely, (a) gamma rays 
and (b) fast charged particles (alpha and beta 
rays). Both types are biologically destructive; 
however, the charged particle 1aciations are not 


2 Neutrons are not emitted by fallout particles. 


very penetrating. The fast beta-rays (electrons) 
can produce burns on unprotected skin, especially 
when fallout material is in contact with the skin; 
but very little shielding, such as is provided by 
average to heavy clothing, is capable of providing 
almost complete protection. The alpha rays are 
not even able to penetrate through the external 
layer of skin. In problems of ingestion of con- 
taminated food, water, and air, however, both 
alpha and beta rays must be considered, because 
when taken internally their full energy is effective 
m causing damage to the internal tissues [16]. 

The gamma rays present the primary shielding 
problem because they are extremely penetrating 
as well as biologically destructive. They may be 
deseribed as high energy X-rays, but for historical 
reasons the term '*y-rays" is frequently used when 
the radiation is produeed by nuclear processes 
and “X-rays?” when it is produced Бу eleetron 
bombardment of an anode. 


8. Characteristics of Fallout Gamma Rays 


Both the quality and the quantity of fallout 
gamma rays have been investigated in many 
research studies [17—23]. 

Gamma radiation may be visualized as a stream 
of individual “energy packets," called **photons." 
The energy content of each packet may be re- 
ferred to as the “photon energy.” Quantitative 
descriptions of fallout gamma rays may utilize 
directly the total photon energy crossing unit area 
of a surface per second; alternately they may 
involve a standard detector whose response is in 
some way deseriptive of the strength of the beam, 
as discussed in Sections 9 and 15. 

Qualitative descriptions are nearly always in 
terms of the “spectrum,” which gives the relative 
importanee of the different gamma ray photon 
energies present. A typical spectrum will be a 
list of numbers or a graph giving, e.g., the fraction 
of the total beam energy m the form of photon 
energies between 1 and 2 Mev, between 2 and 3 
Mev, etc. 

Radioactive material is a mixture of many types 
of atoms, each emitting characteristic photon 
energies. These radioactive atoms originated in 
one of two ways: (1) Some were originally part 
of an atom which fissioned in the nuclear explo- 
sion (mixed fission products); (2) others were part of 
the bomb or were nearby and were made radio- 
active through eapture of à neutron (neutron in- 
duced aetivities). “The mixed fission. produets 
have usually been considered the most important 
radiation source, to the extent that the neutron- 
indueed activities have often been neglected in 
calculations pertaming to fallout. From both 
sources a variety of photon energies 1s generated 
over a range extending up to perhaps 3 Mev. 


9. Exposure to Gamma Rays 


Biologieal effeets depend upon the energy ab- 
sorbed from the gamma rays per unit mass of 


biological material [24]. This absorbed energy is 
usually measured in “rads”? (1 rad=100 ergs/ 
отат). 

The energy absorbed in biological material de- 
pends in turn on the strength of the radiation 
field and the duration of exposure. For protection 
purposes, this is measured in terms of the cumu- 
lative effect of the field on a standard material, 
namely air, in units of roentgeus (r) or milli- 
roentgens (mr). 

Though the roentgen is a measure of exposure 
to a radiation field, whereas the rad is a measure 
of the result of this exposure, it has long been 
common practice to blur this distinction by speak- 
ing of the biological effect of a “dose” of y roent- 
gens, meaning the effect resulting from energy 
absorbed during an exposure totalling у roentgens. 
The International Commission on Radiological 
Units has recommended that the term ‘‘exposure 
dose" be used in this context to prevent confusion 
with the energy absorbed in the biological (or 
other) material.’ 

The strength of a gamma ray field at a given 
time is frequently measured іп terms of the ‘‘dose- 
rate," in units of roentgens per hour (7/hr) or, in 
weak fields, in units of milliroentgens per hour 
(ттт). 

Perhaps Ше best way to gain a feel for Ше mag- 
nitude of different exposure dose rates and total 
doses is by noting the effects upon human beings, 
insofar as this is known or surmised. Specific 
effects depend upon the health of the individual 
the time-duration of exposure, etc., but some gen- 
eralized statements about the effects caused by 
short-duration exposures can be made, asin Table 
9.1 below. 


TABLE 9.1. Probable Effects of Fallout Gamma Radiation 
on Humans ° 


Short Term 
Whole-Body 
Exposure 
(roentgens) 


Probable Effect 


0 to 100 
100 to 200 
200 to 300 
300 to 600 
Over 600 


No obvious effects 

Some sickness 

Sickness and some deaths 

Severe sickness and many deaths 
Few, if any, survivors 





а Note that the effects given all refer to the extent of 
general incapacitation following exposure. Effects on the 
gonads, the eyes, ete., which may be extremely important 
but not immediately obvious, have not been included. 


10. Decrease of Fallout Radioactivity With 
Time 


When fallout material has stopped falling upon 
a given location gamma ray intensities produced 


5 Exposure dose of X- or gamma radiation at a certain place 15 a measure of 
the radiation that is based upon its ability to produce ionization. The unit 
of exposure dose of X- or gamma radiation is the roentgen (7). One roentgen 
is an exposure dose of (X or) gamma radiation such that the associated cor- 
puscular emission per 0.001293 g (ram) of air produces, in air, ions carrying 1 
electrostatic unit of quantity of electricity of either sign [24]. 

4 See footnote 3, above. 


by the fallout can be observed to decrease with 
time. If the fallout is not being spatially displaced 
by winds, this decrease is described roughly by 
the following formula: 

D (ty) = (ъъ) PD (t), (10.1) 
where tı, Ь represent different elapsed times fol- 
lowing the explosion and D(t) is the dose-rate at 
time £ [6, Sections 9.4-9.17]. Equation (10.1) 
can be roughly described in another way: For 
every sevenfold increase in.time there is а tenfold 
decrease in radiation intensity. For example, 
if the radiation dose-rate in a given location is 
1,000 r/hr four hours after the explosion, the 
dose-rate twenty-eight hours after the explosion 
would be about 100 r/hr. | 

Another convenient rule-of-thumb relates radia- 
tion dose-rate (r/hr) to total exposure (г): The 
total exposure dose (r) over a very long time 
interval beginning one hour after the explosion 
is numerically about five times the dose-rate 
(r/hr) at one hour in the same place, assuming 
that no displacement or arrival of the fallout 
material has taken place after one hour. "Thus, 
if the dose-rate one hour after the explosion 
were measured to be 30 r/hr, the total exposure 
at the same location over a very long time interval 
following this measurement would be about 
1507. А more general expression is 

E(tıyt2)= 5t DG) — (65/8) 7], (10.2) 
where Е(%,5) is the total exposure dose accumulated 
between times # and t: following an explosion. 
Note that if &=1 and k~ œo, БЕ о 58 
in agreement with the special rule-of-thumb 
mentioned. These rules do not apply if the 
shielding conditions are changed during the time 
interval. 

The functions of the ratio (5/1) appearing in 
eqs (10.1) and (10.2) are given in figure 10.1. 
It should be remembered that the tume variation 
expressed by the foregoing expressions 18 only ap- 
proximate, that it depends upon the predominance 
of fission product radiation over induced radiation, 
and that it becomes a very poor approximation 
when times after the explosion stretch out into many 
months or years [25]. 


EXAMPLE. Тһе dose-rate of radiation in a 
given place 4 hours after an explosion is measured 
to be 32 r/hr. (a) What will the dose-rate be 
48 hours after the burst; and (b) what will be 
the total accumulated exposure dose during the 
time interval from 4 to 48 hours? 

(а) Taking &=4 hours, t;=48 hours, one ob- 
tains ¢,/f,=12. According to figure 10.1, 
(t/t) 12—(12)-12—0.053. Multiplying this by 
П(4)--32 r/hr, one obtains D(t;)—1.7 r/hr at 48 
hours. ANSWER. 


5 For a recent investigation of this time variation see the report USNRDI- 
TR-425, by Carl Miller (U.S. Naval Radiological Defense Laboratory, 
May 27, 1960). 
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although not shown, the lower portions of the 
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inates the spectrum about a week following the 
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The intense, 


changes with time, and this affects the spectrum 


emitted. 
among gamma rays with photon energies froin 


about 0.5 Mev to 2.5 Mev. 
however, most of the radiation is made up of 


photons below 1 Mev. With further increase in 


el 
energies again become very important. 


in figure 11.1. 
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Figure 10.1. Functions for obtaining the radiation dose rates and cumulative dose at different times after fission. 
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11. Time Variation of the Fallout Spectrum 
The fallout gainma ray spectrum changes with 


(b) To obtain the accumulated exposure dose, 


` one determines from figure 10.1 that [1— (2/4) °] 
ша) ^—0.39. The product (5)Х (4 ћг)х (32 


are contained in the fallout mateiral, and each 


goes through а succession of transformations. 
Thus the very nature of the radiating material 


r/hr) X (0.39) 
time [17-23, 
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Figure 11.1. The relative intensities of different spectral components at several times after fission. The height of each box 
¿s proportional to the (energy content) of gamma rays in the energy interval, and the central line indicates an energy assumed 
for all photons їп the interval for purposes of calculation. The lines do not correspond to spectral energies actually present. 
Тһе total intensity is unity for each spectrum. (Ref. 18, Table 11. Volatile components have been removed from the fission 


products.) 


fission is due to an isotope of Lanthanum, which this monograph correspond to the spectrum exist- 
is the fourth member of a radioactive family which ing 1.12 hours after the fission process. This spec- 


begins with Xenon fission fragment. 


The spectral trum is the one given in figure 11.1 (top curve), 


energy of these gamma rays is 1.6 Mev, and they with no additions for neutron capture components. ` 


are unportant over a time interval 
two weeks. 


of more than This choice results from two considerations: (1) A 
large part of the exposure to radiation is apt to 


Because of the time variations of the gamma occur during the first few hours; and (2) the spec- 
ray spectrum, it is necessary to make some decision | trum one hour after the explosion is representative 
regarding choice of a spectrum to be used as the of other early times in the sense of having about 
basis for shielding calculations. It turns out that Ње same penetrability. 
penetration properties of fallout gamma rays are 


less sensitive to spectral changes than might be 12. Мабп иде of Fallout Gamma Ray In- 


supposed, except for very large penetrations— this 


will be discussed in the next Sections. 


tensities and Exposures 
The data in P 


6 Note that Xenon is a rare gas, so that when this type of fission fragment Precise information about likely fallout dose- 


is initially formed, it will not enter into chemical combination. The Xenon 1 1 - - ; : 
quickly (in a few seconds) disintegrates to produce Cesium, which is very rates 15 not available; but thereis some information 


active chemically and which should become attached to fallout materia. about possible dose-rates produced by large (mega- 


Thus Lanthanum and other daughter products of this chain should contri- 


bute to the fallout spectrum. 


ton) explosions. These vary up to thousands of 


г at early times after the explosion, and could 
result in accumulated exposures up to tens of 
thousands of roentgens for unprotected people, 
even outside blast areas (see fig. 5.1). Further- 
more, there exists strong likelihood of overlapping 
fallout fields, increasing total intensities additively. 
Since a human exposure to a few hundred roent- 
gens is considered very severe, one attempts to 
provide shielding sufficient to reduce radiation in- 
tensities by “protection” factors of the order of 


1000 to insure in practically all cases reduction of 
exposures to an order of magnitude below the dan- 
ger levels. In different types of structures, dif- 
ferent “protection factors" turn out to be feasible. 
For example, in an underground shelter it may be 
relatively easy to incorporate a protection factor 
of 5000, whereas in an above-ground shelter, it 
may be difficult to incorporate more than a factor 
of, say, 100. Figure 5.1 illustrates the value of 
protection factors still lower than 100. 


HI. Description and Measurement of the Radiation 


13. Introductory Remarks 


This part contains a brief outline of terms and 
ideas used in the latter parts of this Monograph. 
Also included is a description of some of the char- 
acteristic physical properties of gamma radiation. 
References are given to other publications con- 
taining a more complete discussion of some of the 
topics. 

14. Physical Properties 


Ganuna ray photons act in many ways as if they 
were tiny “bullets.” They travel in а straight line 
until they “collide” with an atom. When a “соШ- 
sion" takes place the photon may vanish alto- 
gether, or it may ‘‘ricochet” and continue traveling 
in a new direction with lower energy. The ‘‘colli- 
sions" are usually referred to as interactions. 
Those in which the photon disappears altogether 
are called absorptive interactions. ‘The others in 
which the photon “‘ricochets,”’ are called scattering 
interactions [12, Part Ај. 

Аз a photon moves along its trajectory, there 15 
always a possibility that in the next segment of 
the path an interaction of some type may take 
place. Although the nature and location of an 
interaction is a matter of chance, the probability 
of occurrence for each interaction type per unit 
path length traveled depends upon the photon 
energy and the type of material being traversed. 
At very low photon energies an absorptive inter- 
action has such a large probability that a long 
path unimpeded by such an interaction is ex- 
tremely unlikely. 

The average distance between interactions in a 
given type of material is called the “mean free 
path." "The mean free path is proportional to the 
density of the material traversed, and this is the 
main reason why the gamma rays travel about 
700 times farther between interactions in air than 
in water. Mean free paths in air for fallout 
gamma rays are of the order of several hundred 
feet, 1.е., large compared with the dimensions of 
most buildings. This has the consequence that 
most gamma ray interactions within buildings do 
not take place in the air, but in the walls, floors, 
and objects of furniture. 

When gamma rays ''scatter" they continue 
with lower energy in a new direction, so that they 
may undergo one or more additional interactions. 


- 


There is a correlation between the angle of de- 
flection and the energy loss which accompanies a 
scattering interaction. On the basis of this cor- 
relation the following statements can be made: 
Gamma rays deflected through 90° by a scattering 
mteraction cannot have a residual photon energy 
above 0.511 Mev. Similarly, gamma rays which 
have been deflected through 180° by a scattering 
interaction cannot retain an energy above 0.256 
Mev. ‘Thus, gamma rays which have undergone 
large changes in direction are apt to be much 
lower in energy than unscattered gamma rays. 
This is the case even when the direction change is 
the cumulative result of several interactions. 


15. Radiation Detectors 


An object which is affected by radiation to a 
measurable extent can be used as a radiation 
detector [28]. Since the effect of radiation on 
human beings is of primary interest, one might 
wish to discuss all shielding problems in terms 
of tissue damage to the exposed human body. 
This is impractical because the body does not 
make a very suitable detector for experimental or 
theoretical investigations. Instead, detectors are 
used whose response has similarities to the response 
of a huinan body, while being more easily measured 
and interpreted. ‘These detectors tell us about 
the radiation field. "The tissue damage caused by 
exposure to a known radiation field is treated as 
a separate problem, and we do not concern our- 
selves with it further. 

Every detector is characterized by its ‘‘response 
function,” which is the “efficiency” of gamma 
rays of photon energy Е, incident on the detector 
from various directions, in producing the measured 
effect. It is frequently advantageous to have 
an “isotropic” detector, 1.е., one equally efficient 
regardless of the direction of incidence of the 
photons.’ 

The most commonly used radiation detector 
consists of an enclosed air space containing a 
pair of oppositely charged electrodes. When 
exposed to radiation, the gamma rays interact 
with the material surrounding the air space, 
producing high-velocity recoil electrons, some of 


7 Note that the response function depends, in general, on both dircction 
and photon energy and might be denoted by «(Е,д,ф). An isotropic detector 
would correspondingly have a response function e( E), i.e., independent of ê, ф. 


which traverse Ше cavity and ionize the аш 
within. A measurable electron current is pro- 
duced as the ions travel to the electrodes and 
deposit their charge. It is this current, divided 
by the volume of the air space, which constitutes 
the response of the instrument. 

These so-called “cavity ionization chambers" 
are usually constructed in such a way that their 
response is nearly isotropic. But the current 
per emš measured by the instruments is not pro- 
portional to the total energy incident upon the 
detector. Instead it is nearly proportional to the 
energy deposited per gram of detector material, as 
a result of gamma ray interactions occurring 
within the detector. When such detectors have 
cavity walls of a material similar to air in its 
reaction to irradiation, the response can be 
measured in roentgens and referred to as “dose,” 
or more precisely as the “exposure dose," as ex- 
plained in Section 9. 

Reference to the ‘‘detector response" in this 
monograph will always imply an isotropic de- 
tector of this type, 1.e., of air-equivalent material, 
whose response іп roentgens 1s called the “exposure 
dose." 


16. Measurement of the Radiation Flux From 
Particular Directions 


To determine the flux of radiation traversing a 
given region one needs to introduce a suitable 
detector as a probe. For this purpose one of the 
isotropic detectors just described can be used, 
with a response calibrated in roentgens. "The 
detector should be small enough so that all parts 
of it are exposed to the same radiation flux, and 
also so that its presence does not appreciably 
disturb the flux distribution. 

The response measured 1s due to photons which 
traverse the detector material. These photons 
are characterized by the direction of their trajec- 
tory at the time they enter the detector. We 
identify direction by means of polar coordinates 
6,6 which are measured with respect to an 
arbitrary polar axis through the center of the 
detector and a reference half-plane terminating 
along the polar axis. Both “obliquity,” 0, апа 
“azimuth,” ¢, are determined from a line extending 
out from the center of the detector parallel to the 
trajectory but opposite the direction of travel of 
the photon. The angle which this line forms 
with the polar axis is 0, while ó is the angle between 
the reference plane and the plane terminating 
along the polar axis which contains the line. 
(See fig. 16.1.) Any direction of incidence is 
uniquely specified by values for these two direction 
coordinates. 

Note that the coordinate pairs (0,ф) are in 
one-to-one correspondence with points on a sphere 

$ It would seem more natural, perhaps, to identify (0,9) with the direction 
of photon travel rather than its inverse. We would do this if the current of 
radiation were to be our focus of attention. But we will be more directly 
concerned with radiation sources than with radiation fields. IL is therefore 
natural to refer directions to a polar axis extending from the detector out wards 


towards the radiation source, since this is the reference system which is easy 
to use in solid angle calculations of radiation sources. 


of unit radius centered at the detector. We can 
refer to (0,9) either as a direction or as a point 
on such a unit sphere. 

The total detector response D is made up of 
contributions from all (6,9): If D(0,9)sin0d6d 9 
is the contribution to the detector response due 
to gamma rays with incident trajectories between 
0 and 0--40, and between ¢ and ф-аф, the total 
response from all directions 15 expressed by the 
integral 


л 2T 
p= | 5118090 | аф (6,9) 
0 0 


=| е0) | 4006: (16.1) 


We refer to D(6,¢) as the detector dose angular 
(or directional) distribution. 


One reason for introducing D(6,¢) is simply 
that the radiation affecting a detector may be 
analyzed into components coming, say, from a wall 
or a ceiling. If there is negligible air scattering, 
the photon trajectories of importance all intercept 
both detector and wall or ceiling surface of interest. 
We designate this surface by S, and define a func- 
tion 93(д,ф) which is to be unity if a direction (0,9) 
can intercept both surface and detector, but zero 
otherwise. Then if Ds is the detector response to 
radiation from S, 


Ds= | асов) | 4620,90:0,8- (16.2) 


This expression reduces to (16.1) 1f contributions 
from all (0,9) are included. If D(0,9) does not 
depend upon ф, we write (0,9) = (2r) !D(0), and 


9T 
gs (8) = 1/2т ! 4%а«(0,%) 


to obtain 


= |  d(cos 0) D(8)gs (8) (16.3) 


Note that gs(0) 15 the fraction of Ше azimuth con- 
tributing, at obliquity 0 [29, 30]. 

The directions (60,09) which intercept both a given 
surface (S) and a point at the center of a sphere 
of unit radius form a kind of cone with apex at 
the point. (See fig. 16.2.) The aperture of this 
cone 15 measured by the size of the area intercepted 
on the unit sphere; if unit area is intercepted, the 
surface is said to ‘‘subtend a solid angle of one 
steradian at the point." Since the area of the 
unit sphere is 47, the solid angle subtended at any 
point cannot exceed 4r steradians. We shall find 
it convenient to use a different measure of the solid 
angle in this Monograph: Our unit will equal 2r 
steradians. We shall designate the solid angle in 
these units by w, and to avoid confusion will refer 
to о as the “solid angle fraction." 

The solid angle fraction os subtended by a sur- 
face S can be obtained from eq (16.2) by replacing 
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FIGURE 16.1. А detector about to be traversed by a photon. 
The photon path passes through the edge of the detector. 
The polar angles (0,ф) correspond to the line parallel to the FIGURE 16.2. A surface subtending at the center of a detector 
photon trajectory which passes through the reference point a solid angle whose magnitude in steradians is equal to 
at the center of the detector. the shaded area on the unit sphere. 
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FIGURE 16.3. A circular aperture of radius 8 feet, through 
which radiation passes towards a detector placed 4 feet 


from the aperture and on the perpendicular line through 
ats center. 
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D(0,9) by the constant (2т) 1, ог from eq (16.3) 
by replacing (0) with 1. 


as | ( d (cos 0) | aea aso 


= d(cos 0)0« (0). (16.4) 


For a plane surface extending in all directions to 
infinity, with polar axis taken along the perpen- 
dicular to this surface, gs—1 for 0<cos @<1, so 
that both expressions in eq (16.4) give ws=1 for 
this case. Any single plane surface of finite 
extent will subtend an w smaller than 1. 


EXAMPLE. Calculate the detector response D 
at а position opposite the center of a circular 
aperture of radius 3 ft through which radiation 1s 
entering a building, and distant 4 ft from the 
aperture. Assume that the dose angular distri- 
bution from all points of the aperture has the 
constant value D,, that only radiation entering 
through the aperture contributes to the detector 
response (see fig. 16.3), and neglect air scattering. 

We choose as polar axis the direction from the 
detector through the center of the aperture. Then 
the angular distribution of the dose becomes D4 for 
9< 9—0 13/4 and zero for 6 >Ûmax= tan 3/43 
Since бап”! 3/4— cos"! (4/5), we write 


. 1,cos 02 4/5 
95(6,9) = (у сов 94/5. 

Since both gs and D are independent oÍ %, eq 
(16.2) reduces to eq (16.3), and 


1 
јр | d (cos бу D,—1/5 Ds. ANSWER 
4/5 


Note that for constant Па, as in this case, Ds 15 
always given by Daws, 1.е., that ws=1/5 for the 
aperture of this example. 


17. Radiation Sources 


Fallout material has already been referred to 
as the “source” of the fallout gamma rays. We 
often refer to the gamma ray generating material 
as the primary source, because it is a convenience 
at times to treat any interface through which 
radiation emerges as a “source.” If the radiation 
is not actually produced at such an interface, it is 
termed a secondary source. Note that the aperture 
of the preceding example acts as a secondary 
source. 

Fallout material can be spatially distributed 
into almost any geometric configuration con- 
ceivable; but two configurations are of special 
importance because they occur frequently and also 
because many others of real importance can be 
obtained from these two. One of these is called a 
point source and is approximated by a local con- 
centration of fallout material into a volume of 
negligible dimensions. The other is called a plane 
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source and corresponds to a uniform spread of 
fallout materia] in a plane. If the source is as- 
sumed to cover the whole plane, one speaks of the 
source as an infinite plane source; but if the source 
has finite area it 15 called a finite plane source. 

The strength of a radiation source can be de- 
scribed in various ways. It is customary in physics 
to measure the source strength in terms of the 
number of nuclear disintegrations occurring every 
second. А source is said to have a strength of one 
curie if 3.7Х 10? nuclear disintegrations occur per 
second. On the other hand, in radiology one 
speaks of a source which produces an exposure dose 
of one roentgen per hour at a location one meter 
from the source [see, e.g., 31]. Since protection is 
a relative matter, we discuss source strengths in 
this Monograph in terms of the exposure dose they 
produce at fixed positions; and we do not attempt 
to specify source strengths in curies even though 
this is possible. Our practice here is thus akin 
to that just attributed to radiologists. More 
specifically. the real problem in radiation pro- 
tection is to determine how much safer one loca- 
tion is than another. We therefore define a 
standard 'tunprotected" position and measure the 
protection according to reduction of the exposure 
dose below that experienced in the standard 
position. 


18. Protection and Reduction Factors 


We wish to introduce a quantitative measure of 
the protection afforded by structures, by one 
structure relative to another, or one location in a 
structure relative to another. Also, in order to 
gauge the overall effectiveness of protective con- 
struction, a standard “unprotected’’ position is 
needed for comparison. 

Fallout is assumed to cover uniformly all sur- 
faces according to their horizontal projection, as 
already mentioned. From this primary source 
the radiation travels to the detector, wherever it 
may be. A fairly obvious choice for a completely 
unprotected detector location would be above 
ground in a large, open field. This choice has 
been made frequently; but it has two difficulties: 
One is that the detector response depends upon 
the "roughness" of the ground surface; and the 
other is that location of the source at the ground- 
air interface makes it very difficult to calculate 
the detector response accurately even if the spec- 
trum and strength of the primary source are 
completely known. The first difficulty makes it 
hard to “standardize the standard,” so to speak, 
while the second makes it hard to evaluate the 
standard. 

For these reasons we have chosen as a “‘standard 
unprotected position” a detector location three 
feet above a hypothetical source of the same 
character as the fallout on the ground, but at a 
hypothetical smooth, infinite plane “interface” 
with the ground replaced by compressed air of 


the same density? By making the source plane 
ideally smooth we avoid the ground roughness 
difficulty; and by replacing the ground by “сош- 
pressed air” we make an accurate theoretical 
analysis possible in a fairly realistic case and thus 
avoid the second difficulty. To complete the 
description of the “standard unprotected position” 
it is necessary to add that the 3 ft separation is in 
dry air at 76 cm Hg pressure and a temperature 
of 20° centigrade. 

There are at least three additional reasons for 
this choice of reference detector location: (1) It 
gives an extreme, but not unrealistic estimate of 
the dose to which the centroid of the body is 
exposed in an open contaminated field. (2) 
Given the spectrum and strength of the gamma 
rays emitted per unit area of the primary source, 
the reference dose rate can be calculated to about 
2-3 percent accuracy, which is an order of magni- 
tude more accurate than necessary for most 
shielding estimates. (3) It appears easier and 
more natural to consider ground roughness as an 
additional “bonus” shielding from fallout on the 
ground, because the same “bonus” protection 
does not exist against the fallout on top of the 
protective structure. 

We define the “protection factor” P, as Ше 
ratio of the detector response D, in the standard 
unprotected position to the detector response D 
in a protected position, 1.е., 

P=D,/D. (18.1) 


The reciprocal quantity is what is actually deter- 
mined ш structure analysis, and for еазу reference 
it wil be termed the “reduction factor" and 
usually referred to simply as 0/1. Іп this 
monograph По always refers to detector response 
at the “standard unprotected position.” 

Note that the protection factor is not defined in 
terms of a “standard” spectrum, so that it is a 
quantity which will vary, e.g., with time after 
fission, weapon type, etc. This reflects accurately 
the situation because the protection which a struc- 
ture gives will itself vary with circumstances; and 
in all likelihood, the actual variation will cover a 
wider range than given by a variety of different 
estimates. All fission data in this monograph 
correspond to a spectrum calculated for 1.12 hour 
old mixed fission products, as already indicated. 


19. Attenuation of Gamma Rays 


Before proceeding to the discussion of structure 
analysis, we wish to examine briefly what happens 
to gamma rays emitted by a point source em- 
bedded in some material, in order to form а pic- 
ture of their “history” and an appreciation of the 
different factors leading to their attenuation (re- 
duction in strength). [See also 12, Section 5.] 

Very near to the source most of the gamma rays 
are diverging radially away from their point of 
origin. As they diverge from the common source, 

° Note that the density of the ‘‘compresscd air” is immaterial except for 


purposes of visualization, so that the reference configuration is actually an 
,nfinite plane source in an infinite homogeneous medium. (See Section 25.) 
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they also tend to spread apart from each other. 
The effect can be visualized by thinking of con- 
centric spheres surrounding the radiation source, 
with radii designated by r. As gamma rays from 
the source travel outwards they will intercept 
spheres of increasing radii which have an ever 
greater surface area. Even if all the gamma rays 
pass unhindered through the spheres, they will be 
spread ever more thinly the farther out they travel, 
simply because they will be distributed over a 
larger area. This has the consequence that the 
exposure dose measured by a small detector of 
cross sectional area a will decrease with distance 7 
from the source іп proportion to r-2 because Ше 
detector can only respond to a fraction а/4т7? of 
the total flux emerging from the sources. "This is 
the inverse square law of attenuation with distance 
from a source. 

Each gamma ray photon from the source pene- 
trates unhindered through the surrounding mate- 
rial until an interaction occurs, but the interactions 
are chance occurrences. "There is a fixed prob- 
ability determined by the photon energy and type 
of material traversed that in a given segment of 
the path an interaction will occur. "This means 
that in each succeeding layer of material a fraction 
of the photons undergo their first interaction, and 
the number of photons penetrating the layer with- 
out interaction is correspondingly reduced. We 
say that the ‘“‘unscattered beam" undergoes ехро- 
nential attenuation. Successive layers of material 
equally thick reduce the number of unscattered 
photons by the same fraction, assuming the layers 
to be homogeneous and equally dense. Tables of 
attenuation coefficients are given in [32]. 

If the first interaction is a scattering, the gamma 
ray photon may travel further, though with re- 
duced photon energy and changed direction. 
Other interactions follow, separated by path seg- 
ments which tend to become shorter as the photon 
energy becomes smaller. Eventually an absorp- 
tive interaction terminates each photon “history.” 

Many of the gamma ray photons continue to 
increase their distance from the source even after 
a number of scattering interactions. Thus, at 
any fixed distance there will be gamma rays which 
have been scattered once, twice, many times, 
in addition to some which have not been scattered 
at all. Because all “orders of scattering" must 
be considered as giving a contribution to the 
detector response, the exponential attenuation of 
the unscattered gamma rays does not describe 
accurately the law of attenuation due to the 
combined effect of allinteractions. It is customary 
to modify the exponential law (accurate for 
unscattered photons) by a factor called the 
build-up factor to take into account the presence 
of scattered gamma rays. (Tables of build-up 
factors are given in [33].) 

If the primary source is located in the air, 
gamma rays from the source may strike a wall. 
Because gamma rays can have scattering inter- 
actions which completely reverse their direction, 


some of the radiation will be ‘‘backscattered”’ 
from the wall, re-emerging again into the air. 
A great deal of interest attaches to the fraction 
of the gamma ray energy entering the wall which 
comes back out again in this way, 1.е., the “energy 
albedo.” If the different spectral components 
are weighted according to the efficiency of de- 
tection by an air-equivalent detector, an analogous 
"dose albedo" may be determined. We view 


this latter ratio as a type of reduction factor for 
purpose of comparison with the other factors, 
and refer to 1% simply as the albedo [34].! 

Other definitions might be useful, but for the 
most part these four—the inverse square law, 
the exponential factor, the build-up factor, and 
the albedo—are sufficient for discussion purposes. 
Figure 19.1 illustrates the first three factors in a 
particular case. 


IV. Introduction to Structure Shielding Analysis 


20. Basic Approach 


The analysis of structure shielding can best be 
described in connection with a typical, but ele- 
mentary example. Figure 20.1 shows a ‘“‘block- 
house," 1.e., a building with four walls and a roof. 
Fallout rests on the roof and on the ground around 
the structure. We are primarily interested in the 
dose measured at detector position A. 

Radiation contributing to the detector response 
comes from all directions. Because of the low den- 
sity of air, most of the detected radiation moves in 
straight lines to the detector from points of emer- 
gence at the walls. Thus it is natural to speak of 
detector response due to radiation “from the 
walls," “from the roof," “from the window,” ete. 
We expect that radiation ‘from the walls” will 
have little resemblance to radiation ‘‘from the 
roof" in intensity or directional distribution. Ra- 
diation “from the roof" originates from the fallout 
resting on top of the roof, while radiation “from 
the walls" onginates from fallout on the ground, 
possibly a long distance from the structure. These 
marked differences make it reasonable to differ- 
entiate between detector response ‘‘components”’ 
from the different wall surfaces. 

Radiation from a wall ordinarily has entered the 
wall from the outside. Thus when one speaks of 
radiation “from a wall," the reference may ђе 
either to the radiation passing from wall to detec- 
tor or to radiation initially incident upon the ex- 
terior surface of the wall in question. In other 
words, one can classify the detector response com- 
ponents according to interior or exterior surface 
crossed. For our purposes the latteris usually pref- 
erable and can be expressed as follows: If that 
part of the detector response due to radiation en- 
tering the '4"-th distinct wall section from the 
outside is designated D,, the total detector re- 
sponse will be 

1 


where the sum extends over all distinct exterior 
wall surfaces. Such a classification is useful for 
many reasons. ln addition to differences in inci- 
dent radiation characteristies, the attenuation is 
apt to differ widely from one wall section to an- 
other. This has the consequence that the dominant 
source of radiation can be a particular wall or wall 
section. Further, it is usually a straightforward 
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matter to treat backscattering from interior sur- 
faces by a simple corrective factor.! 

We next select a particular wall surface, namely 
that adjacent to detector positions B and О in 
figure 20.1; only radiation entering and leaving 
this wall surface is to be considered for the 
moment. То contribute to the detector response, 
radiation must first pass through the wall and 
then “locate” the detector. The detector re- 
sponse at B divided by the response at position C 
is approximately the reduction factor due to 
interactions within the wall material. If the 
detector is moved away from B towards A, the 
detector response will decrease further. "This 
additional "attenuation" is due to the “inverse 
square” effect described in the preceding Section. 

We distinguish between these two types of 
attenuation and refer to them separately. The 
first, attenuation by interactions with wall ma- 
terial, we shall call barrier reduction due to barrier 
shielding. The second, accompanying a displace- 
ment of the detector farther from the wall, we 
shall call geometry reduction and speak of it as due 
to geometry shielding. The following subsections 
describe these quantities in greater detail and give 
reasoning which leads to their introduction. Here 
it will suffice to say that geometry reduction can 
often be considered as mainly a limited-source 
effect. The wall acts as a secondary source; and 
if this source were constant on a plane extending 
laterally to infinity, there would be no appreciable 
difference in intensity between detector positions 
А and B. The geometry reduction comes about 
primarily because the wall, being finite, subtends 
a smaller solid angle at А tban at В. 


If a reduction В, of the detector response occurs 
between C and В, and а fuither reduction б, 
occurs when the detector is moved from P to A, 
the detector response D, at А may be written 


m ©, В,Рс (20.2) 


where De is the detector response at а reference 


10 Still other types of albedo factor have been defined, using number of 
pbotons in place of energy or dose. Only the “dose albedo” appears in this 
manual, so tbat the use of ''albedo" to refer to dose albedo sbould not be 
confusing. | ) | 

п By classifying according to exterior surfaces one estimates tbat a certain 
fraction of the radiation entering the room from one wall will affect tbe de- 
tector after backscattering from other walls. An estimate of this corrective 
factor is then made with little concern about the detailed bistory of the radia- 
tion or about effects due to radiation entering tbrough the other walls. 
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Figure 20.1. А blockhouse, with fallout on roof and ground. 


position C, which is here shown as just outside 
the wall. 

The dose from a wall section nearly always will 
be written in the form eq (20.2), so that calcula- 
tion of the D; involves determination of appro- 
priate barrier and geometry reduction factors for 
each type of wall surface. 

Though De in this example is the dose just 
outside the exterior wall, the comparison with the 
standard value D, is of most fundamental interest. 
This change requires definition of а modified 
barrier factor 5',— B1Dc/D,, which expresses the 
combined reduction due to attenuation in the wall 
and exterior air. Henceforth barrier factors will 
be of this type, 1.е., ratios of the detector response 
at position adjacent to a barrier to the standard 
response Do. 

The approach described here depends for its 
success on the fact that only —15 percent of the 
radiation energy incident on a surface 15 back- 
scattered. This relatively low albedo insures that 
most of the detector response 15 due to radiation 
which has had contact with only one wall. 


21. Discussion of Barrier and Geometry Re- 
duction Factors 


In the preceding Section a distinction was made 
between two types of reduction in the detector 
response. The reasoning which leads to this 
distinction 1s essentially as follows: 

Most theoretical data derives from the study of 
an especially simple case in which the medium has 
no irregularities (infinite homogeneous medium). 
The source may be concentrated on a plane in the 
medium or localized at a point. Calculations of 
the detector response as a function of distance 
from the source in such an elementary case can 
be done to within a precision of a few percent. 
Because these calculations have internal consist- 
ency and reliability, and also because they 
usually represent а good approximation to the 
reduction in intensity due to actual slab barriers, 
we nearly always find it convenient to consider 
ratios between actual detector response due to 
radiation from a particular wall, say, and detector 
response as estimated by such infinite medium 
calculations. The resulting ratios express nearly 
all the special characteristics of the configuration, 
but they vary slowly as the position or the con- 
figuration is changed. It is possible to estimate 
these ratios with calculations for an infinite 
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medium, since they express effects such as source 
“finiteness.” 12 

Thus, the detector response is expressed as the 
product of two factors: One gives the attenuation 
in the wall material which would occur if the 
source were of a particularly simple type and the 
medium everywhere uniform in density; and the 
other is a ratio which expresses all the other 
features of the actual configuration. For con- 
venience we give these two factors names. The 
name “barrier factor" for the first is а fairly 
obvious choice. The name “geometry factor" is 
chosen for lack of a more accurate term of reason- 
able length. (Note that the term “configuration 
factor" would be more appropriate but is not used 
here because of its length.) 

One more aspect of this development deserves 
comment. In practice, structure shielding anal- 
ysis requires a great variety of standardized 
functions expressible in the form of tables, graphs, 
and possibly empirical formulas. Many of these 
functions are rather naturally classified as geom- 
etry or barrier factors. We make the attempt to 
carry this classification as far as possible, with the 
consequence that standard functions may be 
referred to as “geometry factors" which are not 
obviously related to other functions also referred 
to as geometry factors. The next two sections 
describe some of the independent variables which 
determine these factors both as to type and as to 
numerical value. 


22. Barrier Shielding and Effective Mass 
Thickness (X) 


The shielding provided by a barrier depends 
upon many variables, among which are 

(1) the weight per unit area of the barrier, 

(2) the type of barrier material, 

(3) the gamma ray spectrum, and 

(4) the directional distribution of the radiation 
striking the barrier. 

We discuss items (3) and (4) first: Data will be 
given only for a single fixed spectrum, namely 
that of figure 11.1, (top); variations of the fallout 
gamma spectrum with time after the burst are 
not taken into account. The directional dis- 
tributions of incident gamma rays vary widely, 
so that results for each of the most important 
cases will be identified with a special symbol. 

Items (1) and (2), the type and thickness of the 
barrier, can be treated together because of a 
fortunate circumstance: Nearly all important 
construction materials have atomic numbers low 
enough so that the attenuation is due primarily to 
scattering interactions, which are independent of 
the energy states occupied by the electrons. 
Thus, the attenuation produced by а barrier 
depends almost completely on how many electrons 
it puts in the path of the gamma rays, and this 15 

12 Basically, reduction due to attenuation in barriers is exponential, while 
reduction due to configuration is a power function. We are in fact using 
the same ideas that lead to the introduction of the Build-up Factor concept. 


Estimation of these new ratios requires a knowledge of dose angular 
distributions. 


simply the product of the number of electrons per 
unit volume with the thickness of the barrier. 

Therefore, to measure the effectiveness of bar- 
rier shielding a parameter Х is used which 15 pro- 
portional to A (barrier thickness), and to < Z/A> 
where р is the density and <Z/A> 18 the ratio of 
atomic charge to atomic mass number, averaged 
over the constituent elements of the barrier.” 
(Note that 27 2/47> is proportional to electrons 
per unit volume, the proportionality constant 
being Avogadro’s number.) In particular, we 
define X by the expression 


Ке РА pht (22.1) 


The factor of two is introduced because < Z/A> 
is nearly 0.5 for such important construction 
materials as brick and concrete, so that 2«7Z/A > 
—] for those materials. If the factor 22 Z/A> = 
1 is treated as a dimensionless proportionality 
constant, X can be measured in units lb/ft? (psf), 
since the unit weight is usually given in lbs/ft?, 
thickness in feet, and lb/ft? x ft=lb/ft?. Because 
X 15 nearly always close to the weight per unit 
area, or mass thickness, we call this quantity the 
effective mass thickness. Table 22.1 gives 
2<Z/A> for a number of common materials, 
together with the density of the solid inaterial. 


TABLE 22.1. Values of 2<2/АО> and p 


Material 2< ZA > p, Density in pef 


اپ س 


62.4 
34.0 (average) 
0.076 


. 115 
Concrete : 144 
Soil (depend- 
ing on water 
100 (average) 
480 


710 





a Lead is included strictly for comparison. It should be 
remembered that the data in this monograph do not apply 
to lead because it absorbs rather than scatters the radi- 
ation. See the last paragraph of this section. 


TABLE 22.2. Mass thickness of common building materials 
used in various building components. 


INDIVIDUAL BUILDING MATERIALS 


ltem Thickness | Weight Psf Component 
Inches 
Adober m I 12 116 Wall 
Asbestos board _______ 34 в 1.7 Do. 
Asbestos, corrugated..| __ 4 Roof, wall 
Asbestos shingles ____ 542 1.8 Do. 


Asphaltroofing,3-ply | ........ | 1 Roof 
read 


y. 
Asphalt roofing, 4-ply | ........ 055 Do. 


& gravel. 


13 In particular, <2/А> => (2) fi, where f: is the fraction of the mass 
+ 


due to the ith constituent. | 

и“ Effective mass thickness” and “equivalent thickness of concrete,” a 
term red in the ODM report [7] are different names for essentially the same 
quantity. 


4 


I 
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TABLE 22.2 Mass thickness of common building materials 
used in various building componenis—Continued 


INDIVIDUAL BUILDING MATERIALS 


Item Thickness | Weight Psf Component 
Inches 
Asphaltroofing, 5-ply | ------—- 6.2 Do. 
& gravel. 
Asphalt shingles-.----f _ 274 По. 
Brick i-e 4 39-1 Wall 
ID). EE. 8 79+10 Do 
Е. 12% 115--15 Do. 
Бс _-__ 17 1542-20 Do. 
Clay tileshingles.....] ........ 154-5 Roof 
Clay tile, structural... 4 18 Wall 
Concrete block, hol- 4 30-4 Do. 
OW. 
1079 SEE 8 55+5 Do. 
НОО eee 2 85-10 Do. 
Concrete, reinforced _ - 1 1215 
De — сео 12 144 
Fiber board _..._..-.- 15 0.8 Wall 
Fiber sheathing....... 15 0.9 Do. 
Gy psum block........ 2 9. 52-1. 5 Do. 
Ipo Е  — 4 12. 52-2. 5 Do. 
Gy psum board -_----- 14 221 Wall, ceiling 
Gypsum, sheathing... 15 2.0 Wall 
Marble facing-------- 2 26 Do. 
Plaster, direetly ap- 3⁄4 5 Wall, ceiling 
plied. 
Plaster on fiber lath. _- 1⁄4 5. 0 Do. 
Plaster on gypsum 15 6.0 Го. 
lath. 
Plaster on metal lath__ 34 6.0 Do. 
Plaster on wood lath.. 34 5.0 Го. 
Plaster, solid- -------- 2 20 Wall 
DoE аза 4 30 По. 
16 .0 Do 


d Ceiling 


1 
1 
Plywood, sheathing- - 1.1 Wall, roof 
Slate E ww —-- 3í6 7.3 Roof 
Steel, corrugated, 20 | ______- 2 Roof, wall 
ga. 
Steel panel, 18ga......| ........ 3. 26 Wall, roof 
Stone е 12 130 Wall 
Stone, cast, facing----- 2 24 Do. 
Stucco, metal lath..... 34 9.0 Do. 
Stucco, wood lath.... 34 8.0 Do. 
Terra cotta facing.. ... 1 5.4 Do. 
Wood block, flooring. 34 10 Floor 
Wood finish flooring __ 2542 285 Do. 
оой enee thing ИЕ - 34 2. 5 и rooí 
ood shingles______ Js E 00 
Wood shingles 615 in. | ........ 1.1 Wall 
to weather. 
Wood siding, 8 in. | __ 1.5 Do. 
bevel. 
| Siding О Е 8 2.5 Do. 
rop. 
| 
COMPOSITE MATERIALS 
Бас «Бек. ЭҚЕ oe | Жажа 
ама guau clay 4 60 Wall 
tile. 
Ceramic ЕЕЕ He oe e 9а 
on mortar Бей __ 3⁄4 11 Do. 
zn | 23 Fl 
елшеш oor 
Concrete, reinforced | ———---. | _____ Do. 
ribbed slabs (see 
section on compo- 
nents). 
Marble or terrazzo 4 50 Do. 
on concrete fill. 

Plaster, hollow wall 4 22 Wall 
with stecl studs. - 
Plaster on Suspended E EES 10 Ceiling 

metal lath. 
Wood finish floor...... AI A. o 
on wood sleepers.... TEE CU _ 2 
El weight concrete | 4 | 28 Floor 


TABLE 22.2 Mass thickness of common building materials 
used in various building components—Continued 


DOORS AND GLASS 





Component 


Item Thickness | Weight Psf 


Inches 
34 


4.5 Small buildings 


Door, wood exterior, 
standard 3 ft-0 in. x 
6 ft-8 in. solid core 
flush panel. 

Door, wood interior, 
standard 2 ft-6 in. x 
6 ft-8 in. hollow 
core. 
solid core----------- 

Door, glass, exterior, 
aluminum edge, 
standard 3 ft-0 in. x 
7 ft-0 in. 

Door, glass, exterior 
standard, 3 ft-0 in. 
x 7 ft-0 in. 

Glass, double 
strength. | 

Glass, plate TA | 


Do. 
Large buildings 


Do. 


| Small buildings 


Large buildings 





a Glass. 


Barrier shielding factors are considered func- 
tions of the effective mass thickness, i.e., В-- B(A). 
Use of the concept of effective mass thickness 
makes it possible to perform calculations even 
when barriers are composed of several layers of 
different materials. One merely sums the effec- 
tive mass thickness (X;) of the different layers to 
obtain the total effective mass thickness X= ХХ, 

1 


It should be remembered that the effective mass 
thickness is not pA, the weight per unit area, 
though the two can be numerically equal. Table 
22.2 summarizes effective mass thickness data for 
a variety of common construction materials and 
combinations. This data is taken from the com- 
pilation of [35]. It is important to realize also 
that the mass thickness of a wall, floor, etc., can 
ordinarily be obtained by use of the usual dead 
load tables in engineering handbooks, though for 
very thick walls (several hundred psf), one may 
have to adjust for hollow spaces within the 
construction (see Section 36). 


EXAMPLE: A barrier is composed of wood and 
concrete layers having weights per unit area 5 
lb/ft? and 50 1/2, respectively. Find the 
effective mass thickness for the barrier. 


X=(5 lb/ft?) х1.06+ (50 Ib/ft?)<1.00=55.3 psf. 
ANSWER. 


Figure 22.1 illustrates the dependence of barrier 
attenuation on the mass thickness and atomic 
number. Notice that even steel (atomic number 
26) gives only slightly more protection than con- 
crete of the same effective mass thickness. If 
desired, the ratio between the curves in figure 
22.1 can be apphed as a correction to take into 
account the minor differences between low atomic 
number material; but with composite barriers this 
must be done with care. 

Materials with high atomic numbers, such as 
Pb, give additional protection; but the difference 
is not so great as 1s commonly supposed. To 
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reduce the intensity of 1 Mev radiation by a fac- 
tor of 1,000 requires a concrete barrier having 
Х--305 psf as compared with a Pb barrier hav- 
ing X=190 psf. But the weight per unit area 
is X(2<Z/A >), which for Pb is 190/.791—240 
lb/ft? while for concrete it is 305/1.0—305 1b/ft?. 
Thus nearly % as much Pb (by weight) is required 
às concrete. 


23. Geometry Reduction and Solid Angle 
Fraction 


To identify other important variables on which 
the detector response depends, we consider again 
the blockhouse illustration and in particular the: 
detector response at position À due to the wall 
on the right (fig. 20.1). 

(1) Solid angle fraction: It is clear that the 
detector response will depend on the “apparent” 
size of the wall, as seen from the detector. 
Since “apparent size" is conveniently measured 
in terms of solid angle or solid angle fraction, we 
say that the detector response will depend upon 
the solid angle fraction о subtended at the 
detector by the radiation source. 

(2) Barrier thickness: The thickness of the wall 
(X) will have an effect upon the detector response 
over and above the attenuation which we describe 
by a barrier factor. This is because the direc- 
tional distribution of radiation emerging from the 
wall affects the detector response but is in turn 
affected by the wall thickness. In most cases and 
up to a limit, thicker walls tend to produce direc- 
tional distributions more and more concentrated 
along the perpendicular to the wall surface. 

(3) Wall shape and detector position: The 
dependence on detector position 15 easy to demon- 
strate in the blockhouse illustration by the fact 
that one can find detector positions at the ceiling 
and floor in which the wall subtends the same o, 
but in which the detector response 1s expected to 
differ considerably. Shape effects can also be 
easily demonstrated: À detector 10 ft out from 
the center of the wall will have one response if the 
wall is long and narrow and another if the wall is 
square, though subtending the same solid angle. 

(4) Type of source: The type of primary source, 
as in the case of barrier factors, varies widely 
enough so that we identify each with a special 
symbol. 

(5) Type of wall material: We treat for the 
present only the case of materials of atomic num- 
ber ~30 or below, this range including the most 
important materials commonly employed in con- 
struction. Аз indicated in the preceding Section, 
these low-Z materials have & certain equivalence 
which we utilize. 

In line with the discussion of Sections 20 and 21, 
we consider all effects other than that due to 

15 The detector response could equally well,be said to be a sum of con- 
tributions from each element ofthe source. Following this idea, the detector 
response could be written in terms of the radiation current flowing through 
the differential area elements of the wall. This point of view is equivalent 
to one we adopt, in which the emphasis is on flux directional distributions 
at the detector, provided the radiation current emerging from the wall is 


independent of position on the wall. We make this assumption throughout 
this manuscript. 
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FIGURE 22.1. Attenuation of 1 Mev gamma rays in different materials (plane perpendicular source). 
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150 


lOO 


and H30 are separated by small, nearly constant factors. 


50 


simple attenuation in Ше wall as expressed Бу а 
single “geometry” factor G. Variations due to X 
and w are taken into account by writing 


G=G(X,w). (235) 


То take shape dependence into account ме герге- 
sent actual source shapes as superpositions oÍ 
simpler, and different, source shapes. This makes 
it possible to estimate geometry reduction factors 
for a set of “standard shapes" which can be com- 
bined to approximate most of the important 
shapes encountered in practice (see Part VII). 
Formally, this means that we write 


G(X ,w) = 2 G(X,wi), (23.2) 
where ((Х,о;) is the geometry factor for one of the 
"standard shapes." Note that to insure the 
correct solid angle fraction, the condition 
a= > jc; must be imposed.” 


2 

Although the blockhouse illustration is useful 
for purposes of identifying variables, it does not 
show how one might refer to quite different 
quantities as "geometry factors.” For this рш- 
pose we use a different example: Consider a heavy 
multistory building with fallout on the (fiat) 
roof and detector in the basement. Radiation 
can reach the detector only after passing through 
several heavy floor slabs. To account for barrier 
attenuation and remove it from consideration, we 
may divide the detector response by the expected 
response if the material between source and 
detector were uniformly distributed, and if the 
radiation source extended unchanged in the plane 
of the roof to infinity in all directions (infinite 
plane source). The resulting ratio, which we 
consider a type of geometry factor, would depend 
upon the following things: (1) The extent of the 
roof, which can be expressed in terms of the 
solid angle fraction subtended by the roof at the 
detector; (2) the thickness (X) of the barrier 
between roof and detector; and (8) the shape and 
position of the roof relative to the detector. 
Equations (23.1) and (23.2) are as applicable to 
this case as to the blockhouse wall. The chief 
difference is that w is determined for the fallout 
on the roof, a primary radiation source, in contrast 
to the blockhouse wall example in which w is 
subtended by a secondary radiation source. 


The preceding illustration also makes it clear 
that the estimation of actual barrier and geometry 
factors consists mainly in finding a simpler con- 
figuration, resembling the actual one in important 
respects, which is amenable to a detailed theo- 
retical or experimental analysis. This point is 
elaborated in Parts V and VI. 


24. Brief Summary of the Analysis Procedure 


The preceding discussion may be summarized 
by stating that the dose D measured in a structure 
(such as a blockhouse) will be represented as the 
sum of contributions D, entering through different 
wall sections and ceiling. Each of the D, will be 
the product of a barrier factor B(.X;) and а geom- 
etry factor С(Х,о). Each geometry factor 
consists of one or more terms corresponding to 
“partial surfaces" chosen to represent approxi- 
mately the actual source surface whether primary 
or secondary. Ií the reference detector response 
for the barrier factors is taken to be the standard 
value D, three feet above an infinite, smooth, 
air-air interface having uniform contamination 
equal to that on and around the structure, the 
complete expression for the detector response 
takes the form 


D=D Zi B(X DG (X nw) (24.1) 


Structure configurations exist which call for 
more complicated expressions than this. For 
example, in maze-type geometries the radiation 
may have to turn corners to reach the detector. 
In basements below grade the detector response 
may consist largely of contributions due to gamma 
rays which have scattered downwards from the 
exposed superstructure through intervening floor 
slabs. Such “higher order” contributions give 
terms involving additional factors, but otherwise 
of the same type as the terms in eq (24.1). 

To differentiate between different types of 
data, we will replace “В” and “G” by a variety 
of other symbols, each referring to a single type 
of data, as indicated in the\preceding Sections. 
We discuss next the theoretical methods for 
estimating D's and (5, i.e., the derivation of 
data for applications. 


V. Data 


25. Introductory Comments 


In Part IV a formal approach to shielding 
calculations has been described. This approach 
developed out of a study of simple configurations 
to be described in Part VI; but it has also been 
designed to emphasize the use of quantities which 
can be estimated theoretically. 


16 For an alternative procedure, see Section 49 and references 29, 30, and 41. 
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In Part V we turn to the presentation and 
discussion of data for applications. First we 
give a description of the “basie” data used to 
obtain the functions for applications, so that an 
appreciation of the limitations of the graphs of 
barrier and geometry factors is possible and also 
so that one can get an idea how additional in- 
formation might be obtained. This discussion can 
probably be skipped or skimmed by readers 
interested in a particular application, though the 


“basic? data have some fairly direct applications 
to shielding problems. 

Not all the factors which one might desire can 
be easily obtained theoretically, or even experi- 
mentally. In applications, the graphs of the 
various factors will no doubt be used in cases 
beyond their range of applicability. This makes 
it especially important to give a clear statement 
of their origin and interpretation. Because this 
type of background information is least ambiguous 
when stated mathematically, no attempt is made 
to avoid mathematical terminology. 


26. Basic Data 


To produce varied types of basic data, digital 
computers have proved essential. They are used 
to generate solutions of the integral equations 
describing the transport, diffusion, and energy 
loss of gamma radiation. These “transport 
equations’’ have been solved by “moment 
methods” and by “Monte Carlo methods.” The 
former yield solutions for a source in a medium 
without boundaries, while the latter permit calcu- 
lation of boundary effects [12, Section 7]. 

All tables and graphs presented here have been 
obtained from four types of basic data: (1) Plane 
isotropic source case; (2) point isotropic source 
case; (3) plane oblique source case; and (4) albedo 
results. Тһе first three were obtained using 
moment methods; while the last required Monte 
Carlo calculations. We discuss each type briefly 
here. Appendix A contains an outline of the 
actual moment calculations of the first three data 
types. 

(1) Plane isotropic source case: Consider ап 
infinite, plane source of fallout radiation imbedded 
in an infinite, homogeneous medium. Distance of 
the detector from the source plane is represented 
by d. Alternately, distances may be referred to 
in terms of X, lb/ft? effective mass thickness 
between source and detector. If the medium is 
air (at 20°, 76 cm Hg), the relation between the 
two 15 

4/Х--13.3 ft of air/psf. (26.1) 
At a distance d from the source plane, the radia- 
tion produces a dose angular distribution which 
we refer to as /(d,cos@), where 0 is the obliquity 
measured from the perpendicular from the de- 
tector to the source plane." We fix the scale of 
((d,cos0) so that an isotropic detector which 
registers D, roentgens three feet above the source 
will register dD=D,/(d,cos6)sinédé roentgens due 
to gamma rays striking the detector between 
obliquities 0 and 0+-d0. This implies that? 


И See Section 16, including the footnote. 
. 18 We assume here that the function Да, cos 0) corresponds to penetration 
іп air, Analogous functions corresponding to penetration in other materials 
such as concrete, will also be referred to as, e.g., КХ, cos 8). They will cor- 
respond to the same source strength as in (26.2) (air penetration case), 
but this implies а different value for the integral over all angles. 
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| eos 0) (8 ft, cos 0) — 1. (26.2) 


Computations of these angular distributions 
have usually involved their representation as a 
sum of Legendre polynomials, 


ца, cos 6) => (n+1/2)1,(d)P, (cos 0). (26.3) 


In one type of calculation the first seven or eight 
functions („(4) were determined. An extrapola- 
tion n> was then performed using the fact 
that Да, cos 0) resembles (cos 0) -1 exp (—a/cos 0) for 
cos 0>0 when Ше unscattered component is im- 
portant. In another type of calculation the 
unscattered component was calculated exactly and 
the scattered component determined by sum- 
ming a series of the type eq (26.3) up to seven ог 
eight terms. (See Appendix A.) 

Figure 26.1 illustrates the appearance and 
behavior of /(d,cos0).? The various features of 
these curves have fairly simple interpretations: 
Near the source plane (small d), ¿(d,cos0) is pro- 
portional to (cos 0)^! for angles approaching but 
less than 2/2 (i.e., cos 0—0+). This results in 
large values for directions nearly parallel to the 
source plane, an Important feature of the radiation 
from this type of source. With increasing pene- 
tration the directional distributions become nearly 
isotropic in the forward direction (for 4=150 ft, 
corresponding to about 0.3 mean free paths for a 
photon of average energy). For still greater 
penetrations the distributions become increasingly 
peaked at cos @=1, as the initially dominant oblique 
components are preferentially removed. 

(2) Point isotropic source case: A source con- 
centrated in a region of negligible dimensions 
generates gamma rays which penetrate outwards 
into an infinite, homogeneous medium. At a 
penetration distance d from the source, the dose 
angular distribution is identified by р(4,со5д), 
where @ is measured relative to the axis from 
detector tosource.” Thescale of p(d, cos 6) can be 
fixed conveniently relative to /(d,cos0) by con- 
sidering a small contaminated area A of a plane 
surface located a radial distance 7 from an isotropic 
detector, with A/r?<1. We identify Оур(?,сов0) X 
A/4arr with the dose angular distribution at the 
detector due to this source; but if the whole plane 
surface were similarly contaminated, we would 
want the total detector response to be D, at a 
distance of 3 ft from the plane. To write this 
condition in the form of an integral we take A to 
be an annular sector of width do and length р9ф, as 
shown in figure 26.2. The detector response D, 
at height 3 ft is then given by the integral over 


19 These curves represent actual results for penetration by a 1.12 hr fission 
spectrum in an infinite medium of water, which is so similar to air in regards 
to interaction probabilities that these curves have been scaled to correspond 
to feet of air. 

20 See Section 16, including the footnote. 
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FIGURE 26.2. А small annular sector of a source plane located a distance 
r= N 8? 4- p? ft from a detector 3 ft above the plane. 


obliquities 0 and the total contribution from all 
area elements: 


со or 1 
p= | ate f de |. d (cos 0) Пр(г, cos 0) 


=D 5 | °" p) > (26.4) 


where P(r) = Jd (cos 09) p(r, cos 0). Cancelling D, 
from eq (26.4) we obtain the normalizing condition 
on p(7,cos?) in terms of its integral over all 
obliquities.?! 

The scattered and unscattered components of 
p(r,cos 0) have quite different characteristic fea- 
tures. The unscattered component, p (г, cos 0) is 
concentrated at the angle 0=0, while the scattered 
component, designated p“) (г, cos 0), is distributed 
over all obliquities. As in the case of the plane 
isotropic source, р(г, cos 0) has been calculated by 
using its Legendre polynomial representation.” 
The calculation proceeds more or less as follows: 
A digital computer determines the total strength 
of the unscattered component and also the first 
seven or eight functions p,® (r) by standard 
techniques [12, Part C]. An extrapolation in n 
is then performed, which uses a resemblance of 


‚ 2 As in the case of Ка, cos 6), other functions, corresponding to penetration 
in other materials than air (or water), will be called p( X, cos 0) if the source 
strength is the same as that implied by eq (26.4). 

22 Do not confuse the Legendre polynomial P,(cos 0) with P(r) of eq (26.4). 
Note that P(r) =ро(г). 
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р“ (7,соѕ0) to Ше function (a+cos 0)- 3⁄2. 
Appendix A.) 

Figure 26.3 shows curves of p(r, cos0) for several 
values of r, but rescaled so that the integral of 
pO (r,cos0) over all obliquities is unity. АП 
the distributions are strongly peaked at cos0—1; 
but the truly remarkable feature 1s the similarity 
of trends. About the only change worth noting 
from one value of 7 to the next is in the relative 
strengths of scattered and unscattered compo- 
nents, which is indicated by the separation of the 
curves vertically and (qualitatively) by thelengths 
of the arrows at cosó—1. With increasing pene- 
tration the unscattered copmonent steadily de- 
creases in importance. 

(3) Plane oblique source case: Radiation from 
a plane source goes in all directions; but for some 
purposes it is advantageous to consider only those 
gamma rays which go initially at an obliquity 6 
relative to the normal to the source plane. Тһе 
detector response at different distances from a 
plane source originating gamma rays with the 
fixed obliquity 6) will be designated s(X, созбу), 
with X measuring distance from the source in 
terms of the effective mass thickness (psf). For 
our purpose it is convenient to fix the scale of this 
function so that cos04s9 (X,cos0,) is unity at 
Х--0, where 80 (Х,соѕ0) is the detector re- 
sponse to the unscattered component.” 


(See 


23 Note that Lim cos6os(0 ( X,cos0;) is independent of созбо, for cos@>0. 
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FIGURE 26.4. Gamma ray attenuation curves for monodirectional gamma rays from a plane source, for incident obliquities 
6< 90? relative to the normal to the source plane. (НО, 1.12 hr fission. See also figs. B5 and B6. ) 
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curves of s(X,cos0) for negative соз. Only 
scattered radiation contributes to the detector 
response. The curves are continued to Х--0, 
although strictly speaking the detector response 
changes discontinuously at this point because the 
contribution due to unscattered gamma rays 
begins here. Datafor grazing incidence (cos6,) =0 
is included. The various curves are seen to be 
rather rough. High accuracy is difficult to obtain 
in this type of calculation because a differencing 
of large numbers is involved, and the irregularities 
reflect the inaccuracy of the calculation. 

It 1s by means of the plane oblique data that 
we mast easily make comparisons of the penetra- 
bility and scattering properties of fallout gamma 
rays at different times after fission, together with 
similar comparisons involving (Соб and other 
radioactive materials. Figures 26.6 and 26.7 give 
curves proportional to s(X,1) and s(X,0) for differ- 
ent fission sources, with data also for Соб. In 
figure 26.6, marked differences appear onlv for 
rather large penetrations (X >144 psf). Differ- 
ences in scattering properties are indicated ш fig- 
ure 26.7. The calculations are less reliable, and 
these differences may be exaggerated by the ш- 
accuracy of the data for X near 0. АП curves 
correspond to penetration in water-equivalent 
materials. 

(4) Albedo or backscattering data [34]. The 
gamma ray albedo has been mentioned in Section 
19. If a plane gamma ray source is incident on a 
thick slab, a detector at the interface will observe 
a dose angular distribution D(cos 0). The albedo 
(a) to which we refer 1s defined by 


0 
| d (cos 0) cos 0 D(cos 0) 
|a У а —X (26.6) 
| d (cos 0) cos 0 D (cos 0) 
0 


This ratio depends on the incident source, and for 
inonodirectional sources is a function of the angle 
2 incidence (0) relative to the normal to the slab 
ace. 

Berger and Raso have made machine Monte 
Carlo calculations, for monoenergetic sources, of 
this function [34]. Figure 26.8 gives data for a 
1.12 hr fission spectrum incident on a concrete 
slab, as calculated from their data. 


27. Catalogue of Functions 


The detector response to radiation from an 
arbitrary souice is represented by an expression 
like eqs (16.2) or (16.3), which sums contributions 
from different parts of the source. Source shapes 
can vary enormously, and to cut down the number 
of variables, we do calculations for circular sources, 
illustrated by figure 16.3, with the detector directly 
opposite the center of the circle. Applying eqs 
(16.2) and (16.3) to this simple case, we find that 


Рив) = | d (cos 0) f а.) 
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ОГ 


D(a | М d (сов 6) Ds(6)» (27.1 


where o=1—cos fmax аз discussed in Section 40. 
We approximate both barrier and geometry fac- 
tors by use of integrals of this type with the 
integrand taken as one or another of the functions 
described in Section 26. Interpretation of these 
integrals is given mostly in the next section; here 
we simply list them for intercomparison and 
reference. 

(1) Barrier Factors: The following expressions 
define functions which we use to approximate 
different types of barrier factors: 24 


1, ф= | d (cos 0) l (d, cos Ө), (2/ 
sa= |. d (cos 0) ¿(d, cos 0), (278 

| d (cos 00) s (X, cos 60) 
5 (X)= 2 — (27.4) 

| d (cos 8o) $ (0, cos 0) 
P (d) 35 d (cos 6) p (d, cos 0), (27.5) 
а @= | а (сов 0)p) (d, cos 0), (27.6) 

WX |, d (cos 0) cos 0 s ( X, cos 0) 

Х (1/27) P de l (d, sin 0 cos e). (27.7) 


Most of these quantities have fairly obvious 
designations and interpretations. “Г? stands for 
"layer," “S” stands for '"skyshine," “P” stands 
for "point," and “W” stands for “wall.” Pe 
and p® refer to the scattered component only 
from a point source. 

. (2) Geometry Factors: The ratios to be used 
in estimating geometry factors are listed below: ? 


Г, (Хо) = (1 (Х)-І (=) 
(27.8) 


65 = 1 _ dleos 0) 1 (X, cos 6), 
(27.9) 


24 Notes: P(?(d) is the same as p(? (d) in eq (26.5). Pand РС) should not 
be confused with the Legendre polynomials P,. We consider these defini- 
tions to hold even when the 1, p, and s data corresponds to materials other 
than air (ог H;O). 

25 1{ is assumed that both denominator and numerator in these ratios 
correspond to penetration in the same material, by gamma га уз from the same 
source, 

Note that in the alternative approach mentioned in Section 49 and de- 
scribed in references 29, 30, and 41, it would be logical to use Legendre coeffi- 
cient ratios in a manner closely analogous to the use of the ratios defined here. 
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FIGURE 26.5. Gamma ray attenuation curves 
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The curves are normalized to equal total energy dissipation, integrated over all penetrations. 
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Ly CX o) =) T. d (cos 0) s (LX, cos Ө), 
(27.10) 


DA = [| i _ d (eos 6) (1/27) 


w ае КА, sin 0 cos e), (27.11) 
0 
2 і 
W, (а,ә) -sm |, а d (cos 0) От) 
“| de КА, sin 6 cos e), (27.12) 


S, (dw) E | “4 (сов 0) ( (d, cos 0), 
(27.13) 


P, (d,w) == | | d (cos 0) p (d, cos 0), 
(27.14) 


1 1 
Ре (d o) =P IN | d (cos 0) p? (d, cos 0), 
(27.15) 


PO (do) == | - d(cos б) (1/27) 
T ао р(а, sin 0 cos e). (27.16) 
0 


All these definitions are complete except for eq 
(27.10). The factor v multiplying the right side 
of eq (27.10) 1s chosen so that L,(.X,2)=1, i.e., so 
that the integral over all directions gives the plane 
isotropic source result, L(X). Since L(X), P(X), 
and s(X, cos 0) are connected by the relations 


Cog 


LOO eg | Sor РОС = || КОО гета OY. 


(27.17) 


it turns out that „= (1/2)Р(0).26 

The names assigned to the ratios eqs (27.8) to 
(27.16) are less obvious than in the case of barrier 
factors. In general the letter used, “L”, “Р”, 
etc., is determined by the type of barrier factor 
appearing in the denominator. А more detailed 
description of these quantities is reserved for Ше 
next Section. 


26 This is easily seen for monoenergetic sources, for which 
РО(Х) + P(0)e-2X; s(0( X, cos 0) = (cos 0)-1е-иХ/сов 9, cos 0>0. 


Inserting these functions into eq (27.17), and remembering that the unseat- 
tered component is completely dominant in the limit Х-50, the value for у 
is immediately obtained. But eq (27.17) holds for all .X and for scattered as 
well as unscattercd radiation, since both integrals represent the plane iso- 
tropic dose rigorously. Extension of the argument to include polychromatic 
sources is therefore straightforward. 


622827 0—62——3 
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Notice that in some cases d and in some cases 
X is the penetration variable. The two аге inter- 
changeable, provided the proper conversion is 
used (eq (26.1)). 


28. Description of Graphical Data 


We next take each of the functions given in the 
preceding Section, describe the physical situation 
for which it applies most accurately, present one 
or more graphs, and discuss important features 
exhibited by these graphs. 

u О ГА” int eq (27.2), L(d) ist defined 
as the integral over all obliquities of the quantity 
l(d,cos 0). Therefore L(d) is simply the total de- 
tector response at a distance d (in air) from an 
infinite, plane, isotropic source, divided by the 
total detector response at 3 ft in air from the same 
source. L(X) is the same function except that 
"distance from the source" is expressed in terms 
of the effective mass thickness. Figure 28.1a 1s a 
sketch of the simple configuration for which this 
quantity is directly applicable. "The absorbing 
medium 15 homogeneous, and the detector response 
is isotropic. Figures 28.2a and 28.2b give typical 
curves of this quantity. Both figures show that 
near the source plane the radiation intensity falls 
rapidly with increasing distance from the source. 
This is characteristic of gamma rays having initial 
directions nearly parallel to the source plane (see 
fig. 26.1). After considerable penetration, the 
trend of the curves 15 nearly that for gamma rays 
which leave the source plane perpendicularly (see 
fig. 26.3). 

(2) S(d): The definition of this quantity given 
in eq (27.3) differs from the definition of L(d) in 
regards to the limits of integration. Only obhq- 
uity angles greater than 90? are included. Since 
gamma rays with these obliquities are traveling 
towards the source plane, they have been back- 
scattered. Thus this function gives the total de- 
tector response to skyshine radiation at a distance 
d 1n air from an infinite, plane, isotropic source in 
a homogeneous medium. Figure 28.1b sketches 
the physical situation. Figure 28.3 gives a curve 
of S(d), obtained from the data of figure 26.1. 
The intensity of the backscattered radiation can 
be seen to change rather slowly with distance from 
the source out to perhaps a mean free path of the 
source radiation. As indicated on the graph, 
S (0) =0.098. 

(3) S’(X): This quantity (eq 27.4) is a ratio of 
integrals and is unity for X—0. Both integrals 
sum over negative соб, 1.e., over gamma ray slant 
sources of radiation which travels initially away 
from the detector. The configuration is sketched 
in figure 28.1с, which shows a plane source emitting 
gamma rays isotropically in an upward direction 
into an infinite homogeneous medium. The de- 
tector 1s below the source plane, separated from it 
by au effective mass thickness X. Radiation can 
be backscattered to the detector, giving a response 
measured by the function S'(.X). The normal- 
ization to unity at Х=0 is for convenient use as 


a reduction faetor in the applications. Data are 
given in figure 28.4. The curve bears а consid- 
erable similarity to s(X,0) as given in figure 26.5; 
but this is not surprising, since the grazing inci- 
dence component of thesource contributes strongly 
to this function. 

(4) P(d), P(X): The integral in eq (27.5) in- 
cludes all obliquities, so that P(d) relates to the 
total detector response at a distance d from a point 
isotropic source in an infinite homogeneous me- 
dium (see fig. 28.1d). Results obtained from the 
data of figure 26.2 are given in figures 28.5a and 
28.5b. The scale has been chosen to simplify the 
use of the curve in determining detector response 
due to a small plane area (A) of contamination, as 
discussed in Section 26: 


DID, =P (d) АЈАта?. (28.1) 


(5) The only difference between P(?(X) and 
P(X) 15 omission of unscattered gamma ray con- 
tributions to the integrand. Thus Р®(Х) re- 
lates to the total detector response due to scat- 
tered gamma rays. Graphs of this function and 
also of Р® (Х)-Р(Х)-РФ(Х) are given in 
figure 28.6, as obtained from the data which also 
vielded figure 26.2. 

(6) W(X,d): This is a more complicated quan- 
tity than those preceding. Figure 28.le is a 
sketch of the relevant physical situation. The 
angular distribution of radiation incident on a 
vertical wall is assumed to be the same as if the wall 
were not present. This distribution, /(d,sin@cos¢), 
corresponds to a distance din air from an in- 
finite plane source in a homogeneous medium, 
but with directions referred to a polar axis perpen- 
dicular to the vertical wall (1е., parallel to the 
source plane).*” Assuming the spectrum. incident 
on the wall to be that generated by the primary 
source on the ground, an integral over all oblique 
angles of incidence is made, each weighted accord- 
ing to our information about oblique penetration 
(fig. 26.4). The result approximates the total de- 
tector response at a height d above the source 
plane and with a thickness X of barrier material 
between detector and external barrier surface. 
Errors in this calculation arise from inaccuracies 
in the angular distribution, changes in the incident 
spectrum as the scattered radiation component be- 
comes more important, and incorrect correlation 
of spectrum with direction. In general, this func- 
tion is expected to represent the vertical wall bar- 
rier factor accurately for d small compared with a 
mean free path of the source radiation, i.e., d<1000 
ft. Data are given in figure 28.7. The curves 
have a maximum value of W(0,d) = (1/2)L(d), 
corresponding to radiation from Aa/f an infinite 
plane, namely the half beyond the wall. Notice 
that the bottom curves are not expected to be 
reliable. 


27 Notice that #(4,2050) becomes /(d,cosécoso--sin6sinacoso) relative to a 
second polar axis inclined at an angle = to the original axis, with $ measured 
relative to the plane of the two reference axes. When а-т/2, this reduces to 
l(d,sinócosQ). 
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(7) L.(d,o),L,CX,o): The significance of L;,(X,o) 
is illustrated in figure 28.88: A circular area of 
fallout has its center directly opposite a detector 
and separated by a layer of air of thickness d and 
equivalent mass thickness X. Noting eq (26.4), 
together with the relation с056щщх--1-ш, we 
represent the detector response due to the circular 
area in the form 


X/(0— o) , 
» d 


“or P(X)=L 00—1 (+=). 


(28.2) 


Equation (27.8) is the ratio of the response due to 
the circular area to the response due to fallout on 
the whole plane. Data corresponding to penetra- 
tion in water-equivalent materials are given in fig- 
ure 28.9. While this quantity is easily obtained 
from figures 28.2a and 28.2b by differencing, and 
is included in the data of figure 28.18, the curves 
of figure 28.9 are convenient for comparison 
purposes. 

(8) L,(d,w), Г. (Хо): Theintegralin eq (27.9) is 
over a cone of obliquities of incidence on the detector. 
Thus, this ratio expresses the reduction in detector 
response which occurs if an isotropic detector 
(separated from an infinite plane source of fallout 
gamma rays by a distance d in alr, or a correspond- 
ing barrier thickness X) is replaced by a detector 
responding only to gamma rays incident within a 
particular cone of directions, as indicated in 
figure 28.8b. The response cone has aperture 
тах and subtends solid angle fraction w at the 
detector, with о--1-- созбрах Curves of this 
ratio are given in figure 28.10, as obtained 
from the data of figure 26.1. Notice that these 
curves go to unity only for w—2, 1.е., including 
backscattered contributions for which 1<о<2, 
and amounting to ~15 percent. 

(9) L,(d,o), Ly(X,w): The integral of eq (27.10) 
is over a cone of initial source obliquities. The 
function L,(d,w) therefore expresses the fractional 
reduction in detector response occurring if an 
(infinite plane) isotropic source is suddenly 
allowed to emit radiation only into a limited cone 
of directions about the perpendicular toward the 
detector (see fig. 28.8c). The cone of emission 
has an aperture брах defined by cos64,,—1-— w, 
where o 15 the solid angle fraction of the cone of 
emission, subtended at the detector. Graphs of 
Lo, obtained from the data of figure 26.4, are 
given in figure 28.11. These bear a noteworthy 
resemblance to the curves for Г, and L, As in 
the case of figure 28.10, the curves go to unity only 
as о->2, although figure 28.9 doesn't give data 
for the region of emission in directions away 
from the detector (1.e., o1). 

(10) Equation (27.11) defines W,(d,w) as a ratio 
whose numerator is an integral over а cone of 
obliquities of incidence on the detector. The cone 
axis is parallel to the source plane and the cone 
angle (@max) is related to the solid angle fraction 


subtended by the aperture by the usual relation 
созбрах (1-о), as in figure 28.84. The de- 
nominator of eq (27.11) has the same interpreta- 
tion but with cone angle бр„х=7/2, so that all 
radiation incident on one side of the detector is 
included in the response. Curves of W,(d,w) 
are presented in figure 28.12, as obtained from the 
data of figure 26.1 for penetration in water- 
equivalent materials. 

(11) In some cases one wants to know the 
amount of skyshine radiation passing through a 





. a a Fa. 


(a) 
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(d) 


ЕтсовЕ 28.1. Simple detector-source-medium arrangements: 


N 





vertical aperture. The function W,(d,w), defined 
by eq (27.12), is presented here because of possible 
applications of this type. The definition of W, is 
like the definition of W, except for integration 
limits which include only skyshine contributions. 
Thus W, refers to a half-cone of skyshine radia- 
tion of aperture 64,4, about an axis parallel to the 
source plane, with cos04,,,—1— о, where w is the 
solid angle fraction subtended by the whole cone 
(see fig. 28.13a). The denominator 15 the total 


(b) 


(c) 





(o) Isotropic detector, plone isotropic source; (b) isotropic detector shielded сп the side toward the source, plone isotropic source; (c) isotropic detector, source '' isotropic" 
only in directions pointing oway from detector; (d) isotropic detector, point isotropic source; (e) 1sotropic detector ond ptone isotropic source: The rodiotion field at 
height d above the primory source is token os a new source at o penetration distance X to the right of the detector. Im all cases the calculations ore for infinite homo- 
geneous medio, and alt cases but (d) olso correspond to medio with ptane density variations. 
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FIGURE 28.2а. Plane source of fallout radiation: The 
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c) detector response ratio D/Do as a function of effective mass 
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thickness separating the detector and the source plane. 


FIGURE 28.2b. Plane source of fallout radiation: The (isotropi 
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See also figs. B17 and B18.) 


(H20, 1.12 hr fission. 


X, pst 


FIGURE 28.4. Attenuation curve for radiation backscattered from а plane source isotropic over one hemisphere only, corre- 
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sponding to skyshine radiation incident on the ground. 
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FIGURE 28.5a. Point source of fallout radiation: P(d) as a function of distance d in air between detector and source. (HO, 
1.12 hr fission. See also figs. B19 and B20.) 
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Figure 28.5b. Point source of fallout radiation: P(X) as а function of effective mass thickness between source and detector. 
(Н2О, 1.12 hr fission. See also figs. B21 and B22.) 
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FIGURE 28.6. Unscattered and scattered components of P(X). 
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FIGURE 28.7. Detector response ratio for effective mass thickness X separating the detector from a plane fallout source of 
angular distribution =Í dél(d,sin6cosQ). (H20, 1.12 hr fission. See also figs. B25 and B26.) 
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See also figs. B27 and B28.) 


etector response due to circular plane sources of fallout radiation, as shown in 


(НО, 1.12 hr fission. 


figure 28.8a. 
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FIGURE 28.10. Geometry factor for detector response іо radiation striking 
shown in figure 28.8b. (H20, 1.12 hr fission. See 
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28.11. Geometry factor for response of an isotropic detector adj 
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cone of directions, as sketched in figure 28.8c. 
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FIGURE 28.12. Geometry factor describing detector response to radiation incident in a limited cone of directions about an axis 
mae ç the primary source plane, at height d, as shown in figure 28.89. (H20, 1.12 hr fission. See also figs. B33 
an 34. 
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FIGURE 28.13. Simple detector-source-medium arrangements: 
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(FPO, 1.12 hr fission. 


, at height d above ground, as sketched in figure 28.136. 
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perpendicular to the primary plane 
See also figs. B37 and B38.) 


FIGURE 28.15. Geometry factor describing detector response to skyshine radiation incident in a limited cone 
an axis 
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as sketched in figure 28.18с. 
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5 


the detector in a limited cone of directions about the line from source to detector 
See also figs. B39 and B40.) 


IGURE 28.16. Geometry factor describing detector response to scattered gamma rays from a point source of radiation, striking 
hr fission. 
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10” 


point source, striking the detector 


from source to detector, as sketched in figure 28.18d. 


See also figs. B41 and B42.) 


in a limited cone of directions about an axis perpendicular to the line 


FIGURE 28.17. Geometry factor describing detector response to fallout ganima rays from a 
(НО, 1.12 hr fission. 
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FIGURE 28.18. Contours of constant L(X)L,(X,w). 
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FIGURE 28.19. Contours of constant L(X) L4(X o). 
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skyshine intensity. Data for this function ob- 
tained from the data of figure 26.1 are given in 
figure 28.14. Notice that w=1 corresponds to an 
aperture passing skyshine from one side only, 1.е., 
only half the total skyshine; this is reason for 
the factor 2 in eq (27.12). The variations with d 
were too small to warrant separate curves for 
different d values. 

(12) S,(d,o): This function is defined by an 
integral over а cone of obliquities of incidence 
pointed away from the source plane (see eq 
(27.13)). Thus only skyshine contributes to the 
detector response. The relevant configuration 1s 
sketched in figure 28.13b. "Phe allowed cone has 
aperture angle брах and subtends solid angle frac- 
tion w=1—coS6max. Lhe denominator 1$ the 
response to all backscattered gamma rays incl- 
dent on the detector. Data, given in figure 28.15, 
have been obtained from the data of figure 26.1. 
Note the steepness of the curves near «o— 1, ex- 
pressing the strong component of the gamma ray 
directional distribution in directions parallel to 
the source plane. 

(13) Р.(Фо) and P&?(d,o): All preceding ratios 
are for plane sources. The P, ratios defined by 
eqs (27.14) and (27.15) are like the La ratios, but 
with reference to a point source. Тһе numerator 
is the response of a detector to scattered radiation 
incident within a cone of directions about the 
radial axis from detector to source (see fig. 28.13c). 
The denominator is the total response of an 
isotropic detector to the scattered radiation. 
Data obtained from figure 26.3 is given in figure 
28.16. Note that the curves go to unity only 
for e—2, so that the backscattered contributions 


are not included in figure 28.16. "The following 


expression relates P, and Р: 


Р (Хо) = == [P (X) +P (X) P9 (Ху). 
(28.4) 


This merely expresses total detector response as 
the sum of responses due to scattered gamma rays 
and to unscattered gamma rays. The latter are 
always concentrated along the radial line con- 
necting source and detector. 

(14) Ре (Хо): For completeness we include 
the point source analogue to the W, (plane source) 
function. The function is defined by eq (27.16) 
and illustrated by figure 28.13d. It has in the 
numerator the response of а detector to scattered 
gamma radiation incident within а cone of 
directions perpendicular to the radial line from 
detector to source. The aperture angle бщах 
is related to the solid angle fraction о by 
созвщах-- 1— o. The denominator represents the 
total response of an isotropic detector to scattered 
radiation, while the numerator never represents 
more than half, hence the factor of 2. Data 15 
presented in figure 28.17, as calculated from the 
basic data of figure 26.3. 

(15) А) (Хо), ПХ) Г. (Хә), ГОХ), ә): 
These combinations of barrier and geometry 
factors come up so frequently in applications that 
contour diagrams have been prepared to permit 
direct determination of the product. Figures 
28.18, 28.19, and 28.20 present the diagrams. 
Radial lines have been superimposed and cali- 
brated to permit accurate interpolation. Use of 
French curves facilitates this interpolation. 


VI. Elementary Structure Types 


29. Comments 


A number of elementary barrier arrangements 
have been the subject of special studies. Because 
of their simplicity they permit a fairly complete 
analysis, both experimentally and theoretically; 
and they appear over and over again in various 
combinations in nearly all structures. One might 
say that complicated structures are a composite 
of these simpler types. As an example, the block- 
house described already in Section 20 combines 
vertical walls (Section 33) with a fallout-covered 
shelter (Section 32). 

In this part we discuss these configurations with 
the help of the data and approach already pre- 
sented. Our intention is to describe them sche- 
matically in a manner which leads to numerical 
values for the detector response. At the same 
time we hope to give an appreciation for the nature 
for the approximations involved, so that a feeling 
for size of possible errors can be achieved. 

Because this Part is intended to be illustrative 
rather than exhaustive, we limit the examples to 
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circular source and barrier shapes. The general 
ization to arbitrary detector position relative to 
arbitrary rectangular shapes is given in Part VII. 


30. Density Interface 


The simplest configurations of real importance 
comprise a radiation source at a density interface, 
with the detector placed some distance above in 
air (see fig. 30.1а). Опе may visualize a vast, 
smooth, level field supporting the radiating 
material. 

(1) Air-air, infinite plane source case: In the 
simplest of these cases, the interface is smooth 
and separates two semi-infinite regions, one of 
air at standard pressure and the other of com- 
pressed air of the same density as earth. We 
have already commented that reference to “сош- 
pressed" air has importance only as aid to visual- 
ization, if the source covers the plane uniformly. 
This is because of a scaling principle which 
guarantees that compression of the material into 
plane layers of different density does not alter 
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Figure 30.1. (a) Plane, isotropic fallout source at a simple 
density interface; (b) fallout source on rough ground. 
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FIGURE 30.2. (a) Detector above ап area cleared of fallout; 
(b) detector below a finite source on the ground. 


the radiation flux anywhere, in intensity, spectrum, 
or angular distribution, if “distance from the 
source" is measured in mean free paths, lbs/ft?, 
gm/cm?, effective mass thickness, or other indi- 
cators of the amount of material between source 
and detector. Thus, this case is just the one 
to which calculations of L(X) and L(d) apply 
precisely, and the detector response as a function 
of distance from the source is simply 


DID, L(dy (30.1) 


25 Note that this scaling principle does not hold for non-uniform or finite 
source spatial distributions. 
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HYPOTHETICAL 
E INTERFACE 


FIGURE 30.3. Geometry of a detector above a circular fallout- 
free area, when the fallout elsewhere 28 considered to be 
buried beneath a hypothetical layer of thickness A. 


(2) Air-earth or air-concrete case, infinite plane 
source: If instead of the location at an interface 
between air and compressed air the source is at 
an interface between air and earth, or air and 
concrete, a modification of the detector response 
is expected. This is because the interaction 
probabilities for gamma rays in earth and con- 
crete differ from those of air, for low energy gamma 
rays. This has the consequence that near an 
air-earth interface the flux of low energy photons 
falls a little below tbe flux near the source in an 
air-air case. The effect is small, a few percent, 
and is due to photons which in earth or concrete 
undergo absorptive interactions at higher energies 
than in air. 

This same effect 1s expected for the flux pene- 
trating downward into earth or concrete. By 
and large, earth and concrete can be treated as 
if they were equivalent to water in penetration 
properties, because the differences in interaction 
probabilities give differences in detector response 
which are significant only for very large penetra- 


tions, or for very low energy photons. For our 
purposes it is sufficient to write 
D[Dyz0.9L CX), (30.2) 


within or at the surface of the semi-infinite earth 
or concrete region, if the L(X) data is for pene- 
tration in H;O, as in figure 28.2b. "The factor 
0.9 allows for the lower flux of low energy photons 
in concrete, and should be omitted if the data 
is for penetration in concrete. For an experiment 
of this Туре see [36]. 

(3) Ground roughness: If the interface between 
earth and air is rough, rather than smooth, a 
reduction in intensity is to be expected. The 
reason for this can be appreciated with the aid of 
figure 30.15. That part of the radiation source 
which falls into ‘‘pockets’” cannot contribute as 
much to the radiation component nearly parallel 
to the source plane. The gamma rays emitted 
nearly parallel to the source plane must penetrate 
a substantial amount of dense material before 
emerging into the air, and are correspondingly 
reduced in intensity. 


One expects all earth surfaces to be rough to 
some extent; but this effect will be larger in the 
case of a ploughed field than a carefully smoothed 
baseball diamond. Investigation of ground 
roughness effects is still in a fairly rudimentary 
state. We note, however, that raising the detec- 
tor height above a smooth interface has a similar 
effect of reducing the contribution of the unscat- 
tered "grazing" component through attenuation 
along the elongated path from the point of genera- 
tion. Because of the importance of some type of 
rule for obtaining intensities and angular distribu- 
tions corrected somehow for ground roughness, 
we tentatively adopt the simple procedure, sug- 
gested by Ksanda, Moshkin, Shapiro, and others, 


of writing 
DID; = L(d-4-7), (30.3) 


where т is a constant which is characteristic of the 
local terrain [37]. Values of 7 are not available 
for widely different terrains. At least one case 
has been investigated experimentally, however, 
and for this case it has been calculated that + z40 
ft.22 The experiments were performed in an open 
Nevada field without special treatment either to 
increase or decrease roughness. Equation (30.3) 
has no real justification beyond the comment that 
introduction of the parameter + affects both in- 
tensity and angular distributions the way ground 
roughness affects them. The argument can be 
made to sound slightly more convincing by de- 
scribing the fallout as effectively buried beneath 
a thin layer of earth. 

(4) Finite plane sources: The preceding dis- 
cussions assume that the source covers the plane of 
the interface uniformly. If the region of contam- 
ination does not cover the source plane, new and 
complicating aspects to the problem appear. 
Consider, for example, a circular clearing in a 
fallout field, with the detector above the center of 
the field (fig. 30.2a). Due to ground roughness, 
the detector response will be far from the value 
expected for a smooth source plane. Hardly any 
unscattered radiation will reach the detector 
because it will be blocked on its path along the 
ground by irregularities. 

Another complicating feature of this type of 
problem involves the scattered component of the 
radiation. Gamma rays which travel long dis- 
tances on a path which is always not far above 
ground have a large probability of being deflected 
into the ground. When this happens, distances 
between successive interactions shrink from hun- 
dreds of feet to inches. If the photon is to travel 
any significant distance farther, it must re-emerge 
into the air. .This effect reduces the detector 
response in just the sort of configuration illus- 
trated in figure 30.2a by an amount which is not 
well known but which depends on the radius of 
the clearing. Studies of this “path forcsborten- 
ing" effect have been made, both theoretically 
and experimentally; but results for clearings are 
not available. 


2 C. M. Eisenhauer, private communication. 
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In spite of these difficulties it is possible to 
obtain estimates of a sort of the detector response 
to radiation above a clearing, taking ground 
roughness into account. To accomplish this we 
follow consistently the schematization in which the 
effect of ground roughness is approximated by 
performing a calculation at an increased height 
above ground. This implies that ground rough- 
uess has effects similar to those due to “burying”’ 
the source under a thin layer of earth, as already 
mentioned. In figure 30.3 a detector is pictured 
at a height d above ground. The source is buried 
under a layer of earth which has effective mass 
thickness A, equivalent to a layer of air of thick- 
ness т. The clearing is circular, with the edge of 
the source at a radius d/tan@ out from the point 
below the detector. The distance in air from the 
detector to the edge of the contamination is 
d/cos@, and the additional penetration through the 
covering layer is A/cos@, which is equivalent to a 
distance in air of +/cos0. The penetration from 
the edge of the source to the detector is thus 





equivalent to шт feet of air. 
с050 


tration distance in Ше expression for the detector 


response above a cleared area, we obtain the 
desired approximation: 


DID, L (С) 


1-о 


IÍ we use this репе- 





(30.4) 


where о-- 1--с05д. No allowance is made in this 
expression for path foreshortening. 

One advantage of this sunple approximation is 
the possibility of using it to calculate anisotropic 
detector responses above clearings; but we don’t 
pursue this here. 

The accuracy of our predictions of detector 
response above a clearing is limited by lack of 
knowledge of ground roughness and path fore- 
shortening; but in a related problem which is 
perhaps of greater practical importance these 
complicating effects are not operative. This 
problem is the calculation of the detector response 
in carth (or concrete) below the source plane, when 
the source is of limited extent. Figure 30.2b 
illustrates this casc. Radiation from the source 
traveling initially upward goes such long distances 
through the air that it doesn’t come back to the 
same locality. Thus the detector response will 
be lower than in cases which are similar but 
without the density interface. But this is a small 
effect even when the barrier thickness between 
source and detector is not very great. Itis possible 
to determine the detector response fairly accurately 
from the expression 

P DT OILA LAO), GIU 
where w 1s the solid angle fraction of the source as 
seen from the detector. (The factor .9 18 a rough 
correction when LCX) апа L,(X,w) represent. HO 
penetration data.) 'lhis expression can be seen to 


include a contribution due to radiation just men- 
tioned as being, in fact, absent. Thus it will be 
an overestimate, particularly for Х #0. 


EXAMPLE 1: Calculate the reduction factor 
due to ground roughness, for a detector 3 ft above 
a contaminated area, if the "roughness" parame- 
ter 7 has a value of 40 ft for this type of surface. 

From eq (30.3) and figure 28.2a we obtain 
L(3+40)=0.52. ANSWER. 


EXAMPLE 2: How high above the ground must 
a detector be raised to reduce its response by a 
factor of 1000? 

Applying eq (30.2), but omitting the factor 0.9 
because of the height above the interface, we see 
from figure 28.2a that L(d)—0.001 for d=2950 
ft. ANSWER. 


EXAMPLE 3: Assuming a ground surface such 
that 740 ft, calculate the radius of a clearing 
sufficient to reduce the detector response to 
D,/1000, at a height of 3 ft above ground. 
Applying eq (30.4), using figure 28.2a, we find 


that z (0 )=0.001 for 3 —2940 ft. 





— From 
== 
this we see that 1—w=3/p=43/2940, giving p= 
3(2940)/43 ~200 ft. ANSWER. 

This last should be considered only the crudest 
sort of estimate, in view of the state of knowledge 
of ground roughness effects. 


31. Foxhole Problems 


A typical foxhole configuration is shown in 
figure 31.1. Radiation generated at the ground- 
air interface can enter the aperture even though 
it faces away from the primary source, by back- 
scattering from the air. This backscattered сош- 
ponent of the radiation is also called the “sky- 
shine". Normally the aperture will permit 
skvshine radiation to strike the detector only if it 
travels in an allowed cone of directions as indi- 
cated in figure 31.1 by the dashed line on the right. 

The configuration bears a close correspondence 
to that of figure 28.13b, in which the detector 15 





FIGURE 31.1. Detector on the centerline of a cylindrical 
foxhole. 
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in an infinite medium but is sensitive only to radia- 
tion incident within an allowed cone. We there- 
fore make estimates of detector response in fox- 
holes by means of the functions S(d) and S;,(d,w), 
which represent infinite medium data relevant to 
the case sketched in figure 28.13b. 

It should be noted that in the foxhole configura- 
tion there will be contributions to the detector 
response due to radiation backscattered from the 
bottom and sides of the foxhole. 'This introduces 
a correction factor of magnitude~1.2. Taking 
this into account, our expression for the skyshine 
detector response becomes 

D[D,—1.28(d--vr)S,(d--v,»), (81.1) 
where d refers to height above ground, + is the 
ground roughness coefficient and w is the solid 
angle fraction of the aperture as seen from the 
detector. 

Note that foxhole-type problems can occur at 
great heights as, for example in the case of a sky- 
ну in the roof of a tall building. (See Section 

5. 

One complicating feature of the foxhole сеоте- 
try is the possibility that radiation may go directly 
from fallout near the “Пр” of the hole to the de- 
tector. Strictly speaking, eq (31.1) applies only 
if the radiation source has been removed to such a 
distance from the edge that it cannot penetrate 
through the earth to the detector. 

Attempts have been made to ascertain whether 
direct radiation from fallout near the lip is signifi- 
cant in comparison with the skyshine radiation. 
This will certainly depend upon the size of the fox- 
hole. The contribution “through the Пр” should 
increase as foxhole size decreases, with the ratio of 
lip contribution to total response approaching 
unity as the foxhole size approaches zero. ‘There 
should also be effects due to foxhole shape: An 
elongated foxhole, or one with the edges of the lip 
rounded off should differ from a cylindrical fox- 
hole, or one with a sharply defined hp. It appears 
that for the ideal case of a foxhole large enough to 
hold a man, with vertical walls and a sharp lip, the 
direct radiation penetrating through the lip is com- 
parable in importance with the skyshine.?? 

Calculations of the lip contribution involve inte- 
grals of the type (see fig. 31.1) 


DID,= | "асов д) ИХ (0), cos 6], (31.2) 


which sums the contributions Írom radiation 
emerging into the foxhole at various depths (X) 
below the ground surface, assuming that the direc- 
попа! distribution as а function of depth is not 
affected by the interface at the foxhole wall. Eq 
(31.2) applies only to a detector on the centerline 
of a cylindrical foxhole; but the generalization to 
other cases is fairly obvious and, in fact, great 
variations are not expected (a) if the detector is 


39 C. М. Eisenhauer, private communication . 


moved off-center or (b) if the foxhole is not cylin- 
drical, so long as the solid angle subtended at the 
detector by the aperture is held constant. 

For most purposes the integral of eq (31.2) can 


be approximated by writing apy. dX 


dX 
d(cos 0) 
dX 
ing to the edge of Ше Пр. This leads to an expres- 
sion for D., dose through the lip: 





a constant value correspond- 





and assigning 


_ о(1—«)(2—) ~ SE. 
pp 9028) | ахих), (81.3) 


in which р is the density of the ground, w= 
1—cos@min 1s the solid angle fraction of the aperture 
and d is the depth below the ground surface. 
Values of the integral in eq (31.3) can be obtained 
from figure 32.2. Clifford, Carruthers, and Cun- 
ningham have very recently reported foxhole ex- 
periments [38]. 


EXAMPLE: Calculate the protection factor pro- 
vided by a cylindrical foxhole 3 ft in diameter and 
4 ft deep at a point on the centerline and 1 ft above 
the bottom of the foxhole. 

The solid angle fraction w=1—3/(3?+ (3/2)?)* 
—0.106. We obtain the skyshine contribution 
using this value and Х=0 in eq (31.1). From 
figure 28.3 we obtain S(0)=0.098, and from figure 
28.15 we obtain S,(0,0.106) 220.044. Thus D/ Do 
=1.2S(0)S,(0,0.106)—=0.0052. 

We calculate the Пр contribution assuming that 
the density of earth is p=100 ре. From figure 
32.2 we obtain using cos0—1—0.106—0.894 the 
value 13.1 for the integral in eq (31.3). Thus 


_ (0.106) (0.894) (1.804) 


D./Do= (1003 (13.1) 20.008. 


Thus the two contributions are of the same order 
of magnitude and the total response is given by 


02/1 =0.013. ANSWER. 
32. Shelter Соуегед with Fallout 


Figure 32.1 is a sketch of a typical shelter cov- 
ered with fallout. А different example of a 
similar configuration is the blockhouse illustration 
in figure 20.1, which has fallout on the roof. In 
both cases the radiation penetrates through the 
roof slab almost as if the source were imbedded 
in an infinite medium. The emerging flux and 
directional distribution both can be well repre- 
sented by infinite medium calculation, and we 
inake use of the analogy between figure 32.1 and 
the directional detector illustrated in figure 28.8b. 
Data for the latter case are employed. 

The detector in figure 32.1 is exposed to БасК- 
scattered radiation and to a total spectrum deter- 
mined for concrete rather than for air. As in the 
foxhole configuration one should presumably in- 
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FIGURE 32.1. Detector оп the centerline ој a cylindrical 
fallout shelter. 


clude a backscattering correction (1.15) and a 
correction (0.9) which takes into account the lack 
of low energy photons in concrete. These two 
corrections almost neutralize one another; and 
since our schematization is not expected to be 
accurate to 5 percent, we omit them both and write 


DI D = L(X)L,(X,o), (32.1) 


where X is the effective mass thickness of the 
roof slab and о is the solid angle fraction of the 
roof as subtended at the detector. 

It can be seen іп figure 32.1 that radiation can 
contribute to the detector response by entering 
through the basement walls rather than through 
the ceiling. To take account of this one should 
include a correction very similar to the “foxhole 
Пр” correction discussed in the preceding Section. 
In fact, almost identically the same argument can 
be made, and it leads to an expression more gen- 
eral than eq (31.3), namely 


001—9) 2-4) (Ру су! 1— 
DJD,— 2 јаха), (822) 


where о--1--сов0 and d is distance below the roof 
slab. Values for this integral аге to be found in 
figure 32.2. Fortunately, this correction tends to 
be relatively small if the slab is moderately thick. 

Experiments on structures equivalent to source- 
covered fallout shelters have been performed by 
Clarke, Batter, and Kaplan [39]. 


EXAMPLE 1: Calculate the reduction factor in a 
cylindrical basement 10 ft in diameter, covered by 
а roof slab of concrete 1 ft thick, and with the 
detector positioned on the centerline 9 ft below 
the roof slab. 

For this case, о--1--0/(02--52)%--0.126. This 
value of о, and a value Х-< 144 psf can be inserted 
into eq (32.1). From figures 28.2b and 28.10, we 
then obtain (144) —0.0067 and L,(144, 0.126) = 
0.257, so that D/D,— (0.0067) (0.257) =0.00172. 
ANSWER. 

Tt is of interest to calculate the correction due 
to radiation entering the compartment through 
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in foxholes and the wall contribution in shelters. 


See also figs. B49 and B50.) 
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FIGURE 32.2. Data for calculating 


the walls instead of Ше ceiling. From figure 
32.2 we obtain, using с080--1--0.126--0.874, the 
value (0.69) (0.73) =.5 for the Integral in eq (32.2). 
Thus, 


р, (0.874) (0.126) (1.874) 


(144) (9) (0.5) =0.00008. 


The correction is about 5 percent. 


33. Vertical Walls 


The case of a vertical barrier between a detector 
and a fallout source on the ground is well illus- 
trated by the vertical wall in figure 20.1. Radia- 
tion from the source on the ground penetrates 
upward into the air, so that at different heights 
about the ground there are directional distributions 
similar to those given in figure 26.1 for an infinite 
homogeneous medium. The actual distributions 
are only “similar” because of modifications due 
to ground roughness and the perturbing effect 
of air-wall and air-ground interfaces. 

Radiation from fallout penetrates vertical walls 
with surprising ease. This is because penetration 
into the wall depends on incident obliquity, 1.е., 
the angle of incidence relative to an axis perpen- 
dicular to the wall (and therefore parallel to the 
source plane). The main feature of the radiation 
incident on the wall can be appreciated by making 
a simple approximation. Referring to figure 26.1, 
note a high concentration of radiation with совд 20, 
thus travelling nearly parallel to the source plane. 
suppose all radiation should travel parallel to the 
source plane, with directions uniformly distributed 
otherwise. If obliquity angles @' are measured 
relative to an axis which is itself parallel to the 
source plane, the corresponding distribution, 
l(d,cos0'), would be uniform in 0’. In other words, 


[(d,cosé’)d(cosé’) ocd6', 
so that 


1 
, 

[(d,cos0^) ocn! (33.1) 
This illustrates that near the ground (small d) 
the actual function /(d,cos6") should have a peak 
value for 6’ near zero. Since decreasing 6’ 
corresponds to increasing penetrability, the con- 
sequence is that fallout radiation from the ground 
and incident upon vertical walls is very 
penetrating. 

For our purposes it seems reasonable to neglect 
the perturbation in the radiation field due to 
air-wall and air-ground interfaces, and to take 
ground roughness into account in the manner 
indicated by eq (30.3). Thus we apply to this 
problem the data represented by the function 
W(X,d); and we write for the detector response 
at height d above ground and behind a thickness 
X of wall, 


ШІМ ОХ). (33.2) 
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The factor 0.9 corrects for the deficiency of low 
energy photons when the material adjacent to the 
detector is concrete rather than air, water, or 
wood. It should be omitted if the calculated 
W(X,d+7r) already assumes a concrete wall. 
If backscattered photons should make no contri- 
bution, the factor should probably be closer to 
0.85; but in general backscattering will occur, 
even for a configuration like the blockhouse of 
figure 20.1, in which there is no material im- 
mediately behind the detector (Position B). 
Thus the corrective factor 0.9 is to be preferred. 

Equation (33.2) does not apply to the detector 
response some distance behind a wall of limited 
size. For example, the detector response at B 
in figure 20.1 due to radiation from the adjacent 
wall is given by eq (33.2) but the detector response 
at A due to radiation from the same wall contains 
a further reduction, as discussed in Section 20. 
At present we do not have directional distributions 
for radiation emerging from such a vertical wall, 
so that at best we can only prescribe rough guesses 
for this reduction factor. 

Various methods can be used to construct 
functions for estimating the detector response 
behind a vertical wall. We can, for example, 
interpolate between W,(d,w), which is reasonably 
accurate in the limit of zero wall thickness, and 
P, (о, о), which represents an upper limit to 
the geometry factor in the limit of very thick 
walls.. To accomplish this we assume that 
radiation not scattered in the walls contributes 
according to the thin-wall function, while radiation 
scattered in the walls contributes according to the 
thick-wall function. Assuming the proportion of 
unscattered gamma rays to be given roughly by 
the ratio P (X)/P(X), which we designate 
b( X), we arrive at a function Wa, defined by 


Wai X,d, w) =b(X) W,(d,w) 
4-1.15[(1—6 (X) |P,9? (© w). 


The factor (1.15) in this expression is to normalize 
P,® to unity at w=1, to make it comparable 
with Wa. 

According to figure 28.16, the function PO (d,w) 
varies slowly with d (or X); and we use the 
bottom curve (for d—1000 ft) to represent 
Ра? ( о шо) in calculations of Wa, without ex- 
pecting to introduce errors more serious than 
those inherent in eq (33.3). Іп general, eq 
(33.3) 1s expected to give conservative estimates 
of the geometry reduction factor because P,O 
overestimates the factors expected for thick walls. 

Using eq (33.3) in lieu of a better cstimate of 
the wall geometry factors, we represent the 
detector response in the form 


D[D,—0.9WCX,d 4-7) Wa (Х,4--т,о), 


where d and X are height above the ground and 
wall thickness, respectively, and w is the solid 
angle fraction subtended at the detector by the 
wall, and the factor 0.9 estimates the correction 


(33.3) 


(33.4) 


of Н.О wall data to concrete wall data. Equa- 
tion (33.4) is intended to apply only to walls 
nearly circular (or square), with the detector oppo- 
site the center of such a wall. It may certainly 
be greatly in error if applied indiscriminately to 
eccentric wall shapes and to off-center detector 
locations. (Note that a correction for backscatter- 
ing has been included implicitly.) 


EXAMPLE 1: How thick must a vertical wall 
be to reduce the detector response immediately 
behind the wall to 0.0017, if the front of the 
wall faces a fallout field 6 ft below the detector 
and having a ground roughness coefficient 7 —40 ft? 

Making use of eq (33.2), we examine figure 
28.7 to find that value of X for which W(.X,46 ft) — 
0.001/0.9. This value turns out to be approxi- 
mately X=238 psf, which is just under 20 in. 
of concrete. ANSWER. 


EXAMPLE 2: И the wall of the preceding exain- 
ple is square, 12 ft on a side, and if the detector 
is 12 ft behind the wall, opposite its center, what 
wall thickness would be required to reduce the 
detector response to 0.0012? 

This calculation is considerably more compli- 
cated. 'To begin with, 1t is necessary to calculate 
the solid angle fraction subtended by the wall at 
the detector. 'This can be most easily done by 
the method to be outlined in Section 41; but it 
could also be done by replacing the square with 
à circle having the same area. Using the method 
of Section 41, we note that for squares, e=1, and 
for a distance from the square equal to the length 
of a side, ;—2. From figure 41.2 we then read 
off the value w=.13. Inserting this value, to- 
gether with the value (4+7)=46 into eq (33.4), 
we select several values of X which we expect to 
bracket the desired value, namely Х--144, 192, 
and 240. For these values we obtain from 
figures 28.6, 28.7, 28.12, and 28.16 the following 
information: 


X 

144 
192 
240 


b(X) Wa(46’,.13) Ра (2,13) Н/1(Х,46,13) W(X,46) D/Do 

0. 182 0. 18 0. 46 0. 47 0.0085 0, 0040 
. 152 .18 . 46 . 48 . 0014 
1127 ‚18 . 46 . 49 . 00051 


. 003 
. 00105 


Plotting these values for D/D, on double-log paper 
we find by interpolation that D/D$4—0.001 for 
X +207 psf, corresponding to about 17 in. of 
concrete. ANSWER 


34. Light Superstructure with Shielded 
Basement 


Figure 34.1 is a sketch of a structure which is 
similar to the blockhouse of figure 20.1 except 
that walls and roof are thin and the configuration 
includes a shielded basement. Our primary inter- 
est is in the detector response in the basement ; and 
radiation penetrating into the basement may have 
originated either at the roof or on the surrounding 
contaminated ground. 
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Because the roof is thin it is no obstacle to 
radiation. Therefore the floor slab protecting 
the basement is subjected to a beam of gamma 
rays from the roof; and this radiation is incident 
on the floor slab along lines intercepting some 
part of the roof. Penetration through the floor 
slab proceeds as if there were a source located at 
the slab which could emit gamma rays only in 
directions corresponding to the directions of travel 
of radiation from the roof. Moreover, the cone 
of incident directions doesn’t change significantly 
within a few feet of any point on the floor slab. 
For these reasons we consider the actual configura- 
tion to be in close correspondence with the infinite 
medium arrangement of figure 28.8c, and we esti- 
mate the detector response by means of data from 
figure 28.11. If we include a factor 0.9 to account 
for suppression of low energy components in 
concrete as compared with water or air when the 
L(X) data is for penetration in water, the detector 
response (roof contribution only) is obtained from 


DJD,-0.9L CX) Eo CX, o). (34.1) 


Here w is the solid angle fraction of the roof as 
seen from the opposite point on the other side of 
the barrier from the detector, as shown in figure 
94.1. 

In figure 34.1 it is clear that lowering the 
detector from its position adjacent to the slab 
will not change much the response to radiation 
from the roof because the floor slab will still 
subtend a large solid angle fraction. But it is 
easy to visualize configurations of this general 
type in which the detector response is sensitive 
to the distance from the barrier. (As a simple 
example one need only reduce the length and 
width of the structure pictured in figure 34.1 
while keeping the height constant.) One may 
wish to include another factor in the expressions 
for D/D, to account for this variation: 


О = (АА 9) Па (Х е”), or 
—LGOO4 C о) Ра(Х а“), 


depending on whether w is nearly unity eq (84.28) 
or much smaller eq (34.2b). Here, w’ is the solid 
angle fraction of the protecting floor slab as 
subtended at the detector. The factor 0.9, which 
corrects the data to concrete, has been omitted 
from expressions eq (34.2) because backscattered 
radiation has been omitted altogether through 
the special normalization used for La and Pa; 
this backscattered radiation gives an increase of 
about a factor 1.15, thus slightly more than 
canceling the 0.9 factor. If data for concrete is 
used, a backscattering factor is appropriate. 

The second part of this problem has to do with 
radiation which originates at the ground and pene- 
trates “in and down." This is much more difficult 
to analyze. We proceed by trying to account for 
major factors one by one: If the heavy protecting 
floor slab were absent the problem would reduce 


(34.2a) 
(34.2b) 
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FIGURE 34.1. Fallout on and around a light cylindrical superstructure covering 


a shelter. 
cone of aperture 8. 


to а fox-hole problem and the detector response 
would be given by an expression like eq (31.1). 
The floor slab introduces an extra attenuation, 
and the real problem is to estimate this attenuation 
as a function of slab thickness X. То this end 
we note that radiation from ground contamination 
must be scattered downwards to be incident on 
the floor slab, and that such incident radiation is 
not greatly affected by the clear area covered by 
the superstructure, for structures of average size. 
Therefore, we expect this radiation to resemble 
the skyshine both in intensity and directional 
distribution. This suggests use of the function 
S'(X) to estimate the reduction of the detector 
response below that expected for а foxhole 
geometry with no floor slab (see Section 28, 
paragraph (3), and fig. 28.1с): 


D[Dy— S(d--r)S,(d--7,9)S' (X), (34.3) 


where d, т, and о all have the interpretations of the 
foxhole configuration. 

Use of the function S, in eq (34.3) 1s justifiable 
if the floor slab does not greatly modify the sky- 
shine directional distribution, i.e., if the floor 
slab is thin. If X is substantial and w isn’t small, 
it would be better to omit this factor. 
EXAMPLE 1: Consider а light, cylindrical super- 
structure rising 25 ft above ground level, and 
having a diameter of 35 ft. Underneath this 
superstructure there is a basement of depth 10 ft, 
as shown in figure 34.2. Calculate the protection 
[actor at the center of the basement floor with 
and without a 6 in. first-floor slab (.X—72 psf). 
Assume a ground roughness coefficient 7 —40 ft. 


622827 0О—62—__5 
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The radiation incident on the barrier slab is largely confined to a 
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FIGURE 34.2. Light superstructure geometry of Example 1. 


First we calculate the roof contribution for the 
two cases, using w=1—25/ ¥ (25)?+ (17.5)?=0.181 
and w’=1—10/7(10)?+(17.5)?=0.504, together 
with X—72 for the case with floor slab, and using 
w= 1—35/4 (35)?+ (17.5)?=0.106, w’ =1, апа X=0 
for the case without a floor slab. From figure 
28.20 we obtain L(0)L,(0,0.106)=0.023 and 
L(72)L,(72,0.181) =0.0125. Using eq (34.2b) for 
both cases we must next determine appropriate 
values of P,(X,w): We note that P,(0,1)=1, 
and that according to eq (28.4), P,(72,0.504) = 
(0.166) 7! [0.049 4-0.118(0.71)]—.8. Therefore, the 
roof contribution turns out to be 


No floor slab: D/D o=(0.023)1=0.023, 
біп. floor slab: D/Do= (0.0125) (0.8) — 0.010. 


Next we perform the “foxhole”? type calcula- 
tion of eq (84.3). The values S(40 ft)—0.068, 
S,(40 ft, 0.504) —0.29, and 8 (72) —0.027 give 


No floor slab: D/Do= (0.068) (0.29) =0.020, 


біп. floor slab: D/D.= (0.068) 
(0.29) (0.027) —0.00054. 


Combining the two contributions, we obtain 
No floor slab: Р=р,/р= (0.0234-0.020) -1=23, 


біп. floor slab: P=D,/D=(0.010+0.00054) 7 
=95. ANSWER. 


Note (a) that the floor slab is most effective 
against radiation from the ground, and (b) that in 
the absence of the floor slab approximately half 
the detector response is from the ground 
component. 


35. Apertures 


The analogy between gamma rays and visible 
light may be used to clarify the effect of openings 
upon the detector response within a structure. 
Just as doors and windows permit light to enter 
а room from outside, according to their size aud 
location, so also do doors and windows permit 
radiation to enter when fallout covers the sur- 
rounding ground. In one respect this analogy 
breaks down, however: Radiation from fallout 
can penetrate through opaque walls, so that if thin 
they may be almost as transparent to radiation 
as windows. In this case the size and location 
of windows makes very little difference. 

This illustrates that there is a connection 
between effective mass thickness of a barrier and 
the permissible design of openings in the barrier 
from the point of view of radiation protection. 
When radiation enters a structure mostly through 
windows, the advantages of heavy walls are 
negated; and conversely, if the walls are very 
thin the window design matters little. 

Openings of many different types occur normally 
in construction, aud there doesn’t exist much 
information, either experimental or theoretical, 
about aperture effects. Nevertheless, some gen- 
eral rules can be given for estimating their contri- 
bution to the detector response. In this section 
we discuss briefly openings in walls and roof; and 
we distinguish between two different cases: 
Apertures between source aud detector, and 
apertures not between source and detector. 

Figure 35.1 is a sketch of circular windows in a 
blockhouse wall and ceiling, with the detector 
placed opposite the center of each. In both cases 
the window represents a type of surface which 
differs from the main wall or ceiling in that the 
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effective mass thickness is nearly zero and the 
solid angle fraction subtended at the detector is 
apt to be small. 

In the skylight case the detector will receive a 
contribution due to skyshine in addition to a 
contribution from fallout on the skylight. The 
relative sizes of these two contributions is of some 
interest. For purposes of determining the sky- 
shine contribution the skylight can be viewed as a 
special type of foxhole, so that the detector 
response 1s given essentially by eq (31.1): 

По (4) За (а, с), (85.1) 
where а is the height above the surrounding ground 
and w is the solid angle fraction of the skylight 
subtended at the detector. Equation (35.1) 
applies only for very small o1 zero effective mass 
thickness of window. If there is a covering of 
fallout on the skylight, its contribution plus the 
skyshine contribution can be estimated by viewing 
the configuration as an example of a fallout- 
covered shelter and applying the arguments 
leading to eq (82.1): 

DID, LGOLCCE, o). (85.2 
Here о is again the solid angle fraction subtended 
by the skylight at the detector, and X, the effec- 
tive mass thickness of the window, is very small?! 


The case of a window in a vertical wall bears a 
close relation to the directionally dependent 
detector pictured in figure 28.8d. Radiation at 
the window entrance (fig. 35.1) should have а 
directional distribution much like that above an 
infinite plane source in an infinite medium, if we 
consider only gamma rays traveling towards the 
window. The detector behind the window m 
figure 35.1 “sees” only that part of the radiation 
field at the opening which travels in an “allowed”’ 
cone of directions, very much in the manner of the 
detector in figure 28.8d. We therefore apply the 
data obtained using eq (27.11) to obtain the 
detector response to radiation passing through the 
window, just as in the case of a thin vertical wall 


(eq (33.4)): 
DID,z0.9W (X, d)W,(d, o) 


where d is the height of the window above ground, 
w is the solid angle fraction of the window sub- 
tended at the detector, and X is the window 
thickness, which is presumed to be small or zero. 
Note that it may be desirable to take account of 
ground roughness by increasing the value of d, and 
that the factor 0.9 is to be omitted if the W(X,d) 
data 15 for a concrete wall. | 
Figure 85.2 shows openings in partitions which 
lie between the source and the detector. Radiation 
eutering through these openings can reach the 


(35.3) 


31 Backscattering has been omitted in both eqs 35.1 and 35.2. Strictly 
speaking a factor 1.15 should be incorporated in eq 35.2 and a factor perhaps 
as big as 1.3 in eq 351. The latter is larger because the radiation 15 lower in 
energy. Note, however, that if HO data is used a factor of 0.9 ог so 15 also 
needed. 
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Figure 35.1. Apertures in а blockhouse wall and roof. 


detector without penetrating as much material as 
radiation which passes through the partitions to 
the detector. If the partitions are very light, one 
can probably take the openings mto account by 
assigning an average value to the effective mass 
thickness between source and detector. On the 
other hand, the partitions may be thick, so that 
radiation penetrating through the partitions is 
significantly reduced m comparison with that 
passing through an opening. In this case the 
opening can be considered as giving a separate 
term (D) in eq (20.1), in which term the effective 
mass thickness between source and detector Is 
given a value equal to the combined thickness of 
partitions between source and detector along lines 
through the opening, and the solid angle fraction 
is that subtended at the detector by the opening. 
If more than one opening lies between source and 
detector, additional terms may be required, each 
mvolving effective mass thickness and solid angle 
fraetion determined by overlapping ог non- 
overlapping parts of these openings.” 

Figure 35.3 shows an example of openings 
which are not between source and detector. 
If the partitions are heavy, the presence of these 
openings can be expected to modify the detector 
response significantly. For example, tne floor 
opening permits a substantial increase in the 
strength of the radiation field directly below, 
while the wall opening gives ready access laterally 
to the detector. 

Notice that the radiation must turn the corner 
in order to take advantage of these openings, 
so that a seattering must occur at some position 
which is favored in the sense of being below the 
floor openings and (or) ш front of the wall opening. 
This is actually а “maze” problem, and is quite 
different from the case of openings which permit 
radiation to travel directly from source to detector. 


3? We say hcre that the openings “overlap” when a straight line from source 
to detector can pass through more than onc opening. 
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FIGURE 35.2. Apertures directly between source and detector. 
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FIGURE 35.3. Apertures not directly between source and 
detector. 





Methods for treating this case are essentially 
those used in maze problems, which are discussed 
in Section 38. 


EXAMPLE: A radiation detector 1s at the center 
of a cubical blockhouse of side length 10 ft, and 
with very heavy walls. Small, circular windows 
of equal radius are placed in the center of the 
roof and one wall. Compare the contributions 
to the detector response from skyshine and direct 
radiation through the roof window and from 
radiation through the window in the wall. 
For definiteness assume NO and w=0.05. 
We then obtain the following: 
Figures 
Roof, skvshine: D/D =S (10ft) S, (10ft, 0.05) 
— (0.084) (0.021)=0.0018 28. 3, 28. 15 


mo = (29) La (А, 0.05) 
1 9. 19 


Root, direct: 
X0 2 


—0.01 
Wall: DID,= W (0, 5ft) W, (5ft, 0.05) 
= (0.45) (0.125)=0.056. 28.7, 28. 12 


From this we see that 0.056/(0.056--0.0144+-0.002) 
— 0.78 of the detector response is from the window 
in the wall, and that of the remainder only about 
11 percent is due to skyshine. ANSWER. 

Note that backscattering corrections would 
change these figures somewhat. 


36. Compartmentalization 


From the point of view of radiation shielding, 
the single most prominent characteristic of 
ordinary buildings might be termed compart- 
mentalization. Most of the mass of the building 
which may be expected to act as a radiation 
shield is concentrated in the floors and walls, 
which divide the building mto compartments. 

Consider a multi-story building, as schematized 
in figure 36.1, with fallout on the roof and on the 
surrounding ground. Radiation from the roof 
fallout penetrates the building traveling down- 
wards. As this radiation penetrates deeper into 
the building, it must pass through successive 
layers of material, namely the floors. On the 
other hand, depending on the trajectory, the 
radiation may or may not pass through the walls. 
Thus the floors represent a sort of divided shield, 
while the walls act as a collimator which favors 
some trajectories over others. 

Notice that the same effects are present also 
ш the case of radiation originating on the ground 
around the building. Аз this radiation penetrates 
into the building it must pass through a set of 
walls. But the floors and the other set of walls 
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FiGvRE 36.1. Multistory, compartmentalized structure with 
fallout on roof and дгоипа. 





are so oriented that the radiation may miss them 
altogether if it follows an appropriate trajectory. 
Thus one set of walls forms a divided barrier 
while the other partitions lie more or less parallel 
to the direction of penetration and give advantage 
to radiation traveling directly inwards. 

In addition to the compartmentalization of 
buildings into rooms and corridors, many or even 
most types of wall and floor construction involve 
compartmentalization on a smaller scale, because 
of voids introduced to reduce weight. The most 
obvious examples are hollow concrete blocks, 
hollow tile, and hollow bricks. Another example 
ol slightly different type is that of floor joists, 
which have the effect of collimation оп the down- 
ward-penetrating radiation. Such compartmen- 
talization can cause overestimation of the 
protection afforded by a wall. 

The simplest way to take account of compart- 
mentalization in building materials 15 that of using 
approximations based on a more uniform distribu- 
tion of material. For example, in the case of 
compartmentalized walls one might estimate the 
effective mass thickness by using weight/volume 
ratios, thus assuming the material in the wall to be 
spread uniformly through its volume.  Experi- 
ments tend to confirm that such estimates are 
accurate enough for most purposes, largely because 
the barriers so constructed are usually rather 
light-weight. Тһе effective mass thicknesses 
given in Table 22.2 were obtained from weight- 
volume ratios. 

Approximations in which the material is as- 
sumed to be uniformly distributed can also be 
used in large, multi-room structures, depending on 
partition thicknesses. The method is essentially 
the same as for a single compartmentalized wall: 
One calculates the total mass of material between 
detector and roof, say, and assumes it to be dis- 
tributed uniformly between. The accuracy of 
such a simple approximation depends primarily 
on the partition thicknesses. In general, one 
expects it to be accurate when the partitions do 
not individually reduce the radiation intensity by 
as much as a factor of e~2.7. From figure 28.2b 
we find that this means partitions of less than 
about 40 psf effective mass thickness. It should 
be remembered, however, that such a limit 1s, for 
the present, more of a guess than an accurately 
known figure; aud that 1% depends very much on 
the configuration, as will be indicated more 
clearly in the next Section in the discussion of 
other simple types of approximate calculations 
for compartmentalized structures. 


37. Collimation and Divided Barrier Effects in 
Compartmentalized Structures 


In a more detailed analysis of compartmental- 
ization, it is advantageous to consider collimation 
and divided barrier effects separately, since they 
are quite differeut, even though they depend on 
the existence of particular directions of interest, 
such as “downward”? and “inward.” In this 


Section we make the distinction by giving separate 
consideration to partitions placed “parallel” and 
“perpendicular” to the direction of interest. (We 
include floors when we speak of “‘partitions.’’) 

(1) Partitions perpendicular to the direction of 
penetration: When a barrier Is divided into layers 
through which radiation must penetrate succes- 
sively, the mass distribution can usually be ex- 
pected to affect the detector response. But there 
15 one case m which this is not зо. When layers 
of material and a plane source do not terminate, 
but extend unchanged over the whole plane, as 
indicated in Section 30, the dose measured by a 
detector is independent of the particular manner 
in which the division mto layers takes place, or 
even the existence of a layered structure, so long 
as the total effective mass thickness between 
source and detector remains fixed. 

When the radiation source is finite in extent, 
or non-uniformly distributed, the nature of the 
layering must be taken into account. We con- 
sider three extreme cases, namely (a) the total 
barrier mass confined to a layer adjacent to the 
source; (b) the total barrier niass confined to a 
layer adjacent to the detector, and (c) the barrier 
material spread uniformly between source and de- 
tector. Other situations can arise which are com- 
binations, or which can be expected to lie soine- 
where between the extremes just indicated. For 
example, the mass may be concentrated іп a layer 
half-way between source and detector, or it may 
be divided into two layers, one adjacent to the 
source and one adjacent to the detector. 

When the total barrier mass 15 confined to a 
layer adjacent to the source, one has essentially 
the fallout-covered shelter discussed in Section 52. 
and the analysis utilizes data for L,CX,o). 

Cases in which the total barrier mass is confined 
to a layer adjacent to the detector have already 
been discussed in Section 34. They correspond to 
a radiation source emitting only within a restricted 
cone of directions, so that calculations for this 
configuration utilize data for (Хш). 

Distributed barriers have been discussed in Sec- 
tion 30. The finite source cases of interest utilize 
data for the L.CX,o). 

Very general arguments indicate that the de- 
tector readings in these three cases shouldn't differ 
greatly, a factor of two usually being sufficient to 
encompass the ED of values. "This is essen- 
tially because the barrier and geometry reduction 
factors account for the large purt of the barrier 
effectiveness, and these do not AI ce significantly 
from one mass distribution to another. "There- 
fore, these three extreme cases can be used in 
representing the different cases of layering in 
practical problems. We do not expect a slight 
mismatch of schematization and configur ation to 
result in significant errors. 

Notice that if the barrier is finite, as m figure 
37.1, radiation can escape out the sides, lowering 
the observed dose and giving а slight increase to 
the barrier effectiveness. 

(2) Partitions parallel to the direction of репе- 
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FIGURE 37.1. a on a finite divided barrier. 
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FIGURE 37.2. Falloul on a simple siruciure with parallel 
ратапопз. 


tration: Collimation effects are more difficult to 
estunate than divided barrier effects. At present 
neither experimental nor theoretical studies are 
available. The discussion given here is mainly 
confined to limiting cases and is therefore incom- 
plete even though indicating the nature of some 
of the main effects. 

In the compartmentalized strueture of figure 
36.1, consider penetration by radiation from fall- 
out ор the roof. It is clear that collimation effects 
will depend on wall thickness, floor thickness, wall 
spacing and the distance of the detector from the 
source in terms of the number of compartments. 
Various limiting cases can be considered: 

(1) The walls may be very thick, so that radi- 
ation can hardly penetrate into adjac ent compart- 
ments. In this case the detector sees mainly 
radiation whieh has not penetrated through the 
walls, and the problem reduces to one of tlie 
divided barrier problems with a finite source. 

(2) The floors may be very thick: In this case 
the radiation only penetrates through the floors 
if it travels always in a nearly vertical direction. 
This is because radiation traveling at a substantial 
obliquity has slight ehance of getting through the 


barrier, except when the oblique portion of the 
trajectory occurs near the exit surface of Ше bar- 
rier. The very thickness of the floors provides a 
sort of collimation in this way. It is possible, 
though very unlikely, for this type of collimation 
to be stronger than that provided by the walls, 
so that it doesn’t matter when the walls are 
present. 

(3) The walls may be thin, so that radiation 
readily travels from one “channel”? to the next. 
In this case the collimation effect is slight unless 
the radiation can penetrate to distances large 
compared with the “channel width.” 

The most difficult cases concern walls and floors 
which are neither thick nor thin. It may be 
possible in such cases to justify the approxima- 
tion in which all material between source and 
detector is distributed uniformly throughout the 
intervening space. But this approximation can 
lead to an overestimate of the shielding provided 
by the compartment partitions. A safer estimate 
can be obtained as follows: If the radiation cannot 
travel from source to detector without passing 
through parallel partitions, the total effective 
inass thickness of these parallel partitions 1s added 
to that of the perpendicular partitions; and the 
detector response is determined with this aug- 
mented barrier thickness. 

To apply this rule the primary radiation source 
is divided into components according to the num- 
ber and type of parallel partitions which must be 
penetrated, as shown in figure 37.2. Radiation 
fron. that part of the source directly above the 
detector is not shielded by parallel partitions; 
whereas radiation originating far to one side of 
the detector may have to go through several 
partitions to contribute to the detector response. 
The detector response for the simple case pictured 
in figure 37.2 would thus be estimated by 


D=L(X+OHL,CX +1401) E LOS FL (СА 
+L(X+D)L(X T tos), 


where wı, w, wg are the solid angle fractions sub- 
tended by Sı, Ss, S4 at the detector. 

Notice that even when the radiation might be 
expected to penetrate parallel partitions at a 
slant angle we do not use the siant penetration 
thickness, but rather the normal thickness of these 
partitions in equation (37.1). "There are several 
reasons for this, but the weightiest 1s our desire not 
to over-estimate the protection factor seriously. 
We are treatiug parallel portions as if they were 
part of the divided barrier, and this procedure 
under-estimates their effect somewhat. 

This approach to the parallel partition shielding 
mereases the amount of calculation necessary 
because each component requires а separate 
calculation. 


(5/7215 


38. Mazes 


Figure 38.1 shows an elementary maze. То 
pass from source to detector the radiation must 
follow openings, because the barrier is much too 
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FIGURE 38.1. A simple maze with two 90° bends. The 
solid angle fractions о” and о are subtended at the detector 
by the near and far sides of the bend. 


thick to permit shortcuts." At least two 
changes of direction are necessary for the gamma 
rays contributing to the detector response in this 
example. 

Radiation readily passes down the initial pas- 
sageway leading from the source. The intensity 
of the radiation diminishes steadily with distance 
from the source, maintaining a rough proportion- 
ality with the solid angle fraction which the source 
subtends at the detector. At the end of the 
first passageway, as at positions all along, the 
detector response is partially due to scattered 
ganima rays, though mostly to gamma rays coming 
directly from the source. Once the corner has 
been turned, direct radiation can no longer reach 
the detector in significant amounts, and the 
detector response is due to scattered radiation 
alniost entirely. 

It is possible to examine this scattered radiation 
to deternine where it “turned the corner." 
This leads to the conclusion that the wall surface 
at the corner which can “see” both source and 
detector contributes a substantial fraction of the 
scattered ganima rays which go down the second 
corridor. It follows that maze shielding can be 
enhanced somewhat by renioving or recessing 
these scattering areas so that they по longer have 
access to both source and detector, or by covering 
the walls at the corner with a layer of tin or lead 
sufficiently thick to reduce the albedo considerably. 

Experiments now in progress are leading to a 
further analysis of radiation “turning the corner" 
into contributions from individual wall surfaces, 
and a report on this subject by J. С. Le Doux 


33 This dependence upon solid angle would not hold for an anisotropic 
source. In ап extreme case, for a plane-parallel source, the detector response 
is independent of the length of the first corridor. 


and A. Chilton is to be published: The results 
are by no means simple. Apparently three differ- 
ent components are distinguishable: (1) Radiation 
scattering from ceiling, floor, and the end wall of 
the first corridor; (2) radiation scattering from 
the end wall of the second corridor, 1.e. the side 
wall of the first corridor; and (3) radiation ‘‘cut- 
ting the corner." Wall positions contributing to 
each component can be seen using figure 38.1. 
As the detector is moved down the second corri- 
dor, each component is expected to be reduced in 
strength according to the solid angle fraction sub- 
tended by the scattering surface at the detector. 
This means that components (1) and (3) vary 
inversely with the cube of the distance and com- 
ponent (2) with the square of the distance of the 
detector from the corner. On the other hand, 
components (1) and (2) do not depend on the 
absolute dimensions of the maze, while component 
(3) decreases in inverse proportion to the magni- 
tude of an increase of scale. That is, if the maze 
size is doubled without change in shape, and with 
the detector in the same relative position, the 
third component will contribute only half as much 
to the detector response at the corresponding 
position in the larger maze, but the other two 
components are unaffected. The third component 
is very important for small mazes and ducts, but 
relatively unimportant for large mazes. 

This type of analysis suggests a simple expres- 
sion for represeuting the detector response in 
mazes, namely 


VII. Rectangular 


39. Introduction 


Most of the surfaces of interest in radiation 
shielding are rectangular or representable by 
rectangles. On the other hand, it is much easier 
to generate data for circular surface shapes, with 
the detector opposite the center of the circle; and 
many, even most, calculations can be carried out 
with acceptable accuracy by applying circular- 
shape-data to rectangular-shape-walls. There- 
fore, in this and the following sections we discuss 
methods for doing this. In the last section we 
consider the problem of rectangular-source-data 
more generally. 

The use of circular-source-data to obtain the 
detector response to a rectangular-shape-source 
involves several steps and а number of simple 
types of calculations, which we list below: 

(1) Calculation of the solid angle fraction of 
circular sectors; 

(2) determination of the solid angle fraction 
subtended by any rectangle at a point on the 
central axis of symmetry (“centered rectangle" 
саве;) 

3 Private communications пот J. C. Le Doux to C. M. Eisenhauer. 
Added Note: Very recent publications оп this subject are: C. M. Eisenhauer, 
Scattering of Со-60 Gamina Radiation in Air Ducts, NBS Tech. Note 74 
(1В161575), Oct. 1960 (Office of Technical Services, U.S. Dept, of Commerce), 
С. W. Terrell, А. J. deris, R. О. Lyday, D. Sperber, Radiation Streaming 


in Shelter Entranceways, ARF-1188-12, Oct. 1960 (Annour Research 
Foundation, Chicago, Ш.). 
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TABLE 38.1. Albedo coefficients 41, В, for mazes with a 
long first corridor and large, square cross section. w, is 
the solid angle of the entrance as subtended at the first bend. 





Wo A; B, 
0. 007 0. 056 0. 08 
. 063 .04 ri 
‚ 165 054 Ge 
D D, x | 
Па) и, (Ao, + B (о ;— O) | (38.1) 


in a maze with № bends, where w;,w’; are the solid 
angle fractions subtended by the far and near 
ends, respectively, of the ?'th corner at the (¿+ 1)’st 
corner on the detector; and D, is the detector 
response at the midpoint of the first corner. 
(See fig. 38.1. Note that wo 1s taken to be 0.) 
The two terms represent the 77? and r^? com- 
ponents mentioned in the preceding paragraph, 
and A;, 6; can be considered effective albedos for 
these two components. Experiments with single 
bend mazes, using Co-60 and a long first corridor, 
give estimates of A,, В, as listed in Table 38.1.3? 
Corresponding estimates for later bends are not 
available, but for a safe overestimate the values 
Кү ПЕН) 221 can be used. А gost estiniate 
А . A, 1 . 
for +>1 might be А, ELB, В,= В” which 
relies on the assumption that (А/В)=(А/В). 


Source Shapes 


(3) determination of the solid angle fraction 
subtended by any rectangle at any point, in 
terins of a linear combination of the simpler 
“centered rectangle" cases; 

(4) approximation of a “centered rectangle" by 
circles or circular sectors, so that the approxima- 
tion has the same solid angle fraction and perhaps 
the same “elongation? as the rectangle; 

(5) determination of geometry factors for 
circular sectors. 10 15 not difficult to see how these 
different types of calculations can be combined 
into a fairly accurate procedure for utilizing the 
circular data: We represent the general rectangle 
case first as a combination of ‘‘centered rec- 
tangles” and then, by representing each “centered 
rectangle" as a combination of circular sectors, 
we obtain the general case as a combination of 
circular sectors. Then we calculate the geometry 
factor for each circular sector independently and 
combine these results to obtain a value for the 
rectangle. 


40. Solid Angle Calculations for Circles and 
Annular Sectors 


The solid angle fraction о subtended by a сп- 
cular or annular sector 18 easy to calculate if the 
% Private communication from C. M. Eisenhauer, based on analysis of 


the data of Terrell, et al, of the preceding foctnote and as yet unpublished 
data due to F. X. Rizzo and A. Quadrado, Brookhaven National Laboratory. 
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FIGURE 40.1. Coordinates describing circles and circular 
sectors for calculations of the solid angle fraction subtended 
at detector D. 


reference point is on the axis (see fig. 40.1). This 
is the only case we consider here. 
If pis the radius of a circular area, 2 the distance 


of the reference point from the sufrace plane, 


and cos а= 2/22 + p°, then 
(40.1) 


w=] — Cos a. 

If two concentric circles have radii р and 2; and 

subtend solid angle fractions сл and а», then the 

solid angle fraction of the annulus between the 

circles 1s (for we >) 

w == (02 wW = COS ау — COS б». (40.2) 

Similarly, the solid angle fraction of the 

annular sector between circles of radii p; and р», 
and between azimuths ¢, and 2, is (p^ i) 


(40.3) 


о == 1 (o— w1), 


where m is the fraction of the azimuth occupied 
by the sector, i.e., 


1 
== (оә-еі). (40.4) 


EXAMPLE: Calculate the solid angle fraction 
subtended by a 30° sector of the annulus between 
two concentric disks of radii 5 ft and 7 ft, as sub- 
tended at a point on the axis 10 ft from the plane 
of the disks. 

Since cos a= 10/4/72 --10?— 0.820, and cos a= 
10/4125 =0.894, while e;—49;—30?—7/6, we have 

6 

m= and w=0.08333 X (0.894 —0.820) —0.0062. 
ANSWER. 
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41. Solid Angle Fractions for Rectangular 
Surfaces * 


It is particularly fortunate that the solid angle 
fraction for a rectangular surface can always be 
calculated analytically. The calculation is simpli- 
fied if the conventions of Cartesian geometry are 
observed. The notation in the following para- 
graphs is that of figure 41.1, with the origin of 
coordinates in the plane of the rectangle and 
opposite the reference point (usually the detector 
location). The perpendicular distance from the 
plane of the rectangle to the reference point is 
2; and >, yi @=1, 2, 3, 4) are the coordinates of 
the four corners of the rectangle. The coordinate 
system 15 chosen so that the z and y axes are 
parallel to the sides of the rectangle; and the 
corners are numbered in counter-clockwise suc- 
cession. 

Having specified а coordinate system, we next 
list several important parameters, namelv, 


e= |0;/2;|, 
ті |г/‹ | (41.1) 
== yet. 


The first two of these are referred to as “ессеп- 
tricity ratios," while the last is +1, depending on 
the quadrant in which the corner is located. 
Finally, we define the function 


2 t 
т (е,ђ) === tan 1 —————əƏ (41.2 
(є,2) “< ^ "EMIT ) 
which has the symmetrv property 
1 n 
т(е, 1) = T (> 2), (41.3) 


In terms of these quantities, the following gen- 
eral rules hold precisely: 

(1) Reference point opposite the center of the 
rectangle: For this case ессе and n,=n,. The 
solid angle fraction o subtended bx the rectangle 
15 given by 

ш==т (61,71). (41.4) 
Extensive tabulations of r(e,7) have been made 
and are presented in figures 41.2a and 41.2b. 

(2) Reference point opposite one corner: The 
solid angle fraction w is given by 


1 
> т(в,71). (41. 5) 


This result is easily obtained from зушшешу Бу 
combining the rectangle of interest with three 


36 For a discussion of comparable problems in illumination see reference 40. 
Note, however, that in illumination the integrals are somewhat modified 
because the detectors are usually plane surfaces, having a response function 
which is not isotropic but is proportional to the cosine of the angle of inci- 
dence relative to the normal to the detector surface. 





Figure 41.1. Cartesian coordinate system used to describe а rectangle and а 
reference point (or detector), for calculations of the solid angle fraction. 


other identically shaped rectangles in the other 
quadrants (see fig. 41.3). The value of 7 in eq 
(41.5) pertains to the larger rectangle comprising 
the four parts together. 

(3) Reference point opposite the center of one 
edge: This arrangement is shown in figure 41.4. 
The detector can be considered opposite the corner 
of two similar, adjoining rectangles, each with 
eccentricity ratios e, 71, 50 that 


(41.6) 


1 
“то T (€1,71) - 


(4) Reference point оп а plane ol symmetry: 
An example of this case 18 pictured in figure 41.5. 
The detector can be considered at the corner of 
two similar, adjoining pairs of overlapping rec- 
tangles. The combination vields the result 


к= [Чут Cem) + qa (ез,73) |. (41.7) 


Notice that the terms add when the origin of со- 
ordinatesisinside the rectangle and subtract when 
it 15 outside. 

(5) General case: Figure 41.6 illustrates the im- 
portant fact that a symmetric combination of 
rectangles can always be constructed which рег- 
mits determination of the solid angle fraction sub- 
tended by алу rectangular surface at anv reference 
point whatsoever. One obtains therefrom the 
general expression 


«= |01т (6,7) — Got (€2,02) + дат (ез, тв) — gaT (64,4) | . 
(41.8) 


Each term on the right corresponds to the solid 
angle fraction for one of the rectangles with a 
corner at the point opposite the detector. Тће 
terms all add if the origin of coordinates falls 
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inside the rectangle; while two terms add and two 

subtract if the origin 1s outside the rectangle. 
Since the 7(e;,9;) correspond to solid angle frac- 

tions for "centered" rectangles, they are referred 

to frequently. We shall therefore make use of the 

shorter notation 

(41.9) 


Ti™= T (ет) 


for these “partial solid angle fractions.” 


EXAMPLE. Find the solid angle fraction sub- 
tended by a rectangular roof with dimensions 25 
ft by 35 ft at points (a) 24 ft below the center, 
(b) 24 ft below one corner, (c) 24 ft below the 
middle of one of the short sides, (d) 24 ft below 
the roof and 10 ft out from one of the short sides, 
on the plane of symmetry, and (e) 24 ft below the 
roof and half-way between the center and one 
corner. All five positions are shown in figure 41.7. 

Calculation of the solid angle fractions can pro- 
ceed in two, sometimes three ways. For many, or 
most, cases 1t is sufficient to determine the param- 
eters е, 7; and then refer to figure 41.2. Thus, 
for the point designated (a) one has «= (25/2)/ 
185/27 7 ИВ m 424 )/(85/2) 1.371 Опе then 
locates the point having these coordinates on fig- 
ure 41.2, places alongside this point a French curve 
which approximately matches the neighboring con- 
tours, and reads from a nearby radial line the 


value 0.173. ANSWER. 


By comparison, for point (b) one first can de- 
termine the values «—]|25/(—35)1—0.7143, m= 
|(24)/(— 35) =0.6857, and then calculate numeri- 
cally the result о=0.25 (2/т) бап”! (0.7402) = 0.1014. 
ANSWER. 


Notice that in all these calculations the co- 
ordinate system 15 oriented in the same way. For 
pomum E mne е 105/35 —0:357, m —24/35 = 
0.686, one obtains from figure 41.2 the value о= 
0.124. ANSWER. 
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This function is the solid angle fraction subtended by “centered rectangles.” 
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FIGURE 41.2а. Contours of constant r(e,n),e <1. This function is the solid angle fraction subtended by “centered rectangles.’ — 
Continued. 
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This function is the solid angle fraction subtended by “centered rectangles.’ 
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FIGURE 41.4. Reference point (or detector) opposite the center of one 
edge of a rectangle. 






To calculate the solid angle fraction for point 
(d), one first lists 


6 =12.5/(85—10)=0.5 
&—|(—12.5)/(—10)|— 1.25 
т--24/(35--10)--0.96; 41--1 
т--|24/(-10)|--2.4; 4ҙ--1. 
Figure 41.2 then provides the results т(0.5, 0.96)= 
0.216, т(1.25, 2.4)=0.111, so that о=0.5(0.216-- 
0.111)=0.164. More precise numerical evalua- 
tion gives the result o=0.1659. ANSWER. 
Finally, for point (e), one has 
«= 103/4) (25)/(3/4)(—35)|=0.714 
e=|(—1/4)(25)/(8/4) (—35)|=0.238 
FiGuRE 41.5. Reference nl delector lan 
E UI nee 2) ыы в--|(-1/4)(25)/(1/4)(35)|--0.714 
74 





«= |(3/4)(25)/ (1/4) (35)| 2.143 

m= |24/(3/4) (—35)|— 0.914; a= — 1 
m=}; 92=1 
ng=|24/(1/4)(35)|=2.743; q5— — 1 


74— n3; q= 1 


Numerical evaluation of the т(е;,7;), followed by 
combination as indicated in eq (41.8), yields the 
result «=0.1524. А trial calculation using figure 
41.2 as a shortcut gave the less precise value 0.153. 


ANSWER. 
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FIGURE 41.6. Reference point (or delector) opposite arbitrary 
localions: (а) inside and (b) outside а rectangle. 
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FIGURE 41.7. Deleclor positions for the Example, relative 
lo a 23 fl by 35 fl reclangular roof. 
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42. Approximation of Rectangles by Circular 
Sectors 


We next cousider simple methods for represent- 
ing rectangles as combinations of annular and 
circular sectors. Ц this can be accomplished for 
“centered rectangles” the extension to the general 
case is straightforward. Of course, it is in prin- 
ciple possible to make such representations with 
any accuracy desired, by using many small 
annular sectors. But this amounts to a form of 
numerical integration, which we seek to avoid. 
To be acceptable, a representation must use a 
minimum number of sectors. 

Three methods wil be discussed, each illus- 
trated by one of the sketches on the left in figure 
42.1. In the first, the elongation of the rectangle 
is neglected altogether. In the second, the elonga- 
tion 15 accounted for Бу the use of two circular 
sectors. In the third, mainly for use with rec- 
tangles far off center, annular sectors are used. 

(1) Approximation by circles: In some cases, 
eccentricity effects can be shown to be weak, so 
that two surfaces can be considered equivalent if 
they subtend the same solid angle fraction. One 
can then calculate the solid angle fraction т (e,n) of 
a “centered rectangle,” and replace the rectangle 
with a circle subtending w=r(e,n) (see fig. 42.1а). 
In the general case of an off-center rectangle, each 
of the four “centered rectangles” is replaced in 
this way by a circle, and the corresponding ap- 
proximation has the form 


: (42.1) 





| | 
== [32 
Те! 


in analogy with eq (41.8), with the o; representing 
circles. Тһе combination can be viewed as а pair 
of annular sectors whose combined solid angle 
fraction is equal to that subtended by the off- 
center rectangle. 

(2) Approximation by two circular sectors: 
This method attempts to obtain greater accuracy 
through an attempt to take eccentricity effects 
into account. In figure 42.1b the rectangle 15 
indicated by а solid line. The dashed lines indi- 
cate two squares, one having 15 sides equal to the 
length of the rectangle, while the other has its 
sides equal to the width of the rectangle. The 
circular arcs represent parts of circles which, if 
complete, would subtend solid angle fractions 
equal to values for the squares. "The representa- 
Поп consists of the complete smaller circle pius 
enough of the annular ring betwcen the circles to 
give the combination a total solid angle fraction 
equal to that of the rectangle. 

If the solid angle fraction of the rectangle 15 
t(e,n), the solid angle fractions subtended by the 
{wo squares are 7(1,5), and r(1,n/e). [tis easy to 
show that the approximation takes the form 


w= (o + (1—a)w’’, (42.2) 





cos 0 
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FIGURE 42.1. Approximation of а “centered rectangle’ by circles and annual sectors: 


(a) Neglecting elongation (Method 1); (b) taking elongation into account in a simple two-term approximation (Method 2); (c) using an annular sector to represent an 
off-center case (Method 8). "The graphs on the right show qualitatively the circular sector approrimations (dashed lines) to the g: (8) functions (solid lines) which correspond 


to the sketches on the left. 
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where w is the solid angle fraction of Ше (“сеп- 
tered”) rectangle, w’ and w’’ represent circles 
having solid angle fractions 7(1,7) and 7(1,7/e), 
respectively, and 


обе) т луд. 
т(1,7) —r(1,n/e) 


Extension to the general case of off-center rec- 
tangle is accomplished by the use of eq (42.1), 
with each of the о; referring to a sector-represented 
“centered rectangle." 

(3) Approximation by annular sectors: This 
method is for use when the four terms in the 
general case cancel strongly against one another, 
as when the rectangle is far to one side of the origin 
of coordinates. 

Concentric circles passing through the farthest 
and nearest corners of the rectangle subtend 


(42.3) 


I——— } and 
Ут? eH 


respectively. А sector 


solid angle fractions о’ 


ДАР" 1 == 73 

па te? 3-1 
of the annulus between these circles can represent Ше 
rectangle if it has the same solid angle fraction. 
Our representation thus takes the form 
(42.4) 


w=mM(w’—w’’), 
where w’, w’’ represent the two circles, with solid 
angle fractions as indicated, and 
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Vn +e? +1 Vm +e +1 

Since this last method is applied to off-center 


rectangles primarily, there is no reason to consider 
combinations of such representations. 


(42.5) 


43. Calculation of Geometry Factors for 
Rectangles 


As a last step in the calculation of geometry 
factors for rectangles, we replace win the preceding 
section by G(X,w), a geometry factor. We can 
do this because Ше G's combine linearly just like 
the ws. Thus, (42.1) takes the form 


GO» СО СХ, (431) 


where Ше term on Ше left is Ше approximate 
geometry factor for the rectangle, and the terms 
on the right are geometry factors for circles sub- 
tending solid angle fractions w;.*” 

37 Note that eq (43.1) is exact when the G( .X,o;) correspond to the centered 


rectangles determined by the four corners of the off-center rectangle, rather 
than to circular sector approximations to these centered rectangles. 


622827 O—62——6 


Similarly, in method (2) utilizing circular 
sectors, eq (42.2) becomes 


G(X w) =aG(X,w’) + (1—a) а Хе”), 
(43.2) 


where Ше value of а is given by eq (42.3). Сош- 
bination of the four results of the type eq (43.2) 
into an expression of the form eq (48.1) gives the 
circular sector approximation for the general case. 

Finally, the annular sector approximation (3) 
for a geometry factor is obtained by translating 
eq (42.4) into 

G(X о) =m[G(X,w’) — G(X,w’’)}, (438.3) 
where m is given by eq (42.5). 

These methods can be applied to any of the ge- 
ometry factors represented generically by G(X,w), 
with greater accuracy in some cases than in others. 
No attempt has been made to compare system- 
atically the accuracy of the different approxima- 
tion methods in different types of calculations. 


EXAMPLE 1: А rectangular underground shelter 
40 ft x 20 ft is covered with a concrete slab of 
effective mass thickness Х= 150 psf. Fallout on 
the slab emits gamma radiation of intensity sufh- 
clent to give a response D, to a detector 3 ft above 
an idealized plane similarly contaminated. Cal- 
culate the detector response 10 ft below the slab, 
half-way between the center and one corner, using 
the circular approximation. Figure 43.1 illus- 
trates the configuration. 

To determine the detector response one first 
calculates the parameters e;, nı, and т; for each of 
the corners of the roof slab, using figure 41.2. 


1 є; Ni т, qi 
1 30/15 10/15 0. 58 1 
2 130/(—5) | 10/(—5) 1278 
3 |(-10)/(-5) 10/(— 5) . 205 1 
4 |(— 10) /15| 10/15 . 40 = 


Using figure 28.10 we next calculate L, for circles 
having the same solid angle fraction: 


Lo Qi L,(150,w;) 
1 0. 58 08723 
2 . 278 . 49 
3 . 205 . 39 
4 . 40 Ж. 








Figure 43.1. Detector below the roof slab of ап underground 
shelter, as in Examples 1, 2, and 3. 


Finally, we make the combination Гог the general 
case: 


D=D,L(150) (1/4)} >; (-)954,1,(150,е) 


= D,(.0057) (.555) =0.0032D). ANSWER. 


The circular approximation has fairly obvious 
limitations. It cannot be applied accurately for 
length/width ratios which differ by more than 
perhaps a factor of three from unity, unless the 
directional distribution involved is either nearly 
isotropic or nearly monodirectional. In the 
example just worked, the term 71.(150,») 15 
expected to be rather inaccurate because во 15 
6; but this term contributes only about 10 percent 
of the final result. 

Clearly this method cau be expected to become 
more inaccurate when the terms cancel rather 
than add, 1.e., if the detector is to one side of the 
roof slab. Also, values obtained will represent 
an overestimate or underestimate depending on 
whether radiation contributes more at angles 
near perpendicularity or at angles nearly grazing 
the source plane. Thus the value obtained in 
this example is an overestimate because radiatiou 
emerges from the roof slab with greater intensity 
at angles near the perpendicular. 


EXAMPLE 2: Re-calculate the detector response 
for the preceding example, using a circular sector 
approximation. 

Using the е, по and т; values already deter- 
mined, the following table is easy to calculate: 








í а? ы?” а 1-a 1,4(150,щ») ћа(150,а7 ") 
1 0.48 0.71 0. 569 0. 431 0. 68 0.78 
2 . 128 CI . 743 . 297 . 26 . 78 
3 . 128 . 336 . 630 . 370 . 26 . 55 
4 . 336 . 48 . 556 . 444 ‚55 ‚ 68 


From these values, using eq (43.2), we obtain 


0 о; Га (150,42) 
1 0. 58 0. 728 
2 . 278 . 394 
3 . 205 . 367 
4 . 40 . 608 


The final combination then yields 
4 

D=D,L(150) (1/4) 2 (—)**!g,L,(150,o;) 
:=1 - 


—D,(0.0057)(0.523)—0.0030D,. ANSWER. 


The value is to be compared with 0.0032D,, 
which was obtained by the circular method. The 
agreement 15 surprisinglv good, considering that 
Ше length/width ratio of the roof slab is 2. The 
largest part of the disagreement appears to come 
from the term corresponding to the second corner, 
with its length/width ratio of 6. The ratio of 
values given by the two approximations for this 
term 1s about 0.81. 
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EXAMPLE 3: Consider the rectangular roof slab 
of the preceding examples as two adjacent squares, 
and calculate the contribution to the detector 
response from the square farthest from the 
detector. "The detector is still considered to be 
10 ft below the roof, halfway between center and 
corner. 

The values of the parameters e;, по and т; turn 
out to have the same magnitudes as before; the 
ond changes are in the q;, which affect the value 
of m: 





2 е: ni т: EC 
1 30/15 10/15 0. 58 

2 P 10/(—5) . 278 — | 
3 10/(—5) 10/(— 5) . 205 -4 
4 10/15 10/15 . 40 1 


The sign changes lead to increased cancellation, 
so that the resulting solid angle fraction is reduced 
to 0.063. 

The values of w and o” сап now be calculated: 


10/15 


5 
(720% 10g . 7 
JG) to 2 w =l 


Inserting these values into eq (42.5), we obtain 
т=.063(21/8) =.165. 
Ме next obtain the values 


, 


=| 


1„(150,ю)--.78, Г. 150 о p n 


from figure 28.10, and make the final combination 
mdicated by eq (43.3), 


D= D,L(150) (0.165) (0.78— 0.54) =0.00023D). 
ANSWER. 


This result is much smaller than the value 
obtained iu the preceding examples partly because 
the solid angle fraction is much smaller and partly 
because the radiation tends to go downward 
rather than sideways after emergence from the 
roof slab. 


44. Further Discussion of Geometry Factors 


It 1s possible to calculate geometry factors while 
avoiding the use of approximations like those de- 
scribed in the preceding Sections. То accomplish this, 
the circular source data obtained from expressions 
typified by eq (27.1) must be replaced bv calcu- 
lations using more general integral forms of types 
eqs (16.2) and (16.3), with 0,(0,е) and g,(0) cor- 
responding to rectangles rather than circles.*5 

The “centered rectangle" case is the only one 
that need be considered in this connection be- 
cause off-center cases are always obtainable as 
linear combinations. Because е and п uniquely 


38 For an alternative and different approach to the calculation of geometry 
factors see Section 49 and references 29, 30, and 41. 


specify a centered rectangle configuration, it is in 
principle possible to make extensive tabulations 
of geometry factors G(X,¢,n), using these variables 
in place of о. Presumably @(X,e,n) can be pre- 
sented, for each X, by means of contour diagrams 
similar to figure 41.2. In the future this may 
prove to be a desirable way of presenting geometry 
factor data, particularly for radiation directional 
distributions which depend on an azimuthal angle 
e as well as a polar angle 0, as in the case of radia- 
tion emerging from a vertical wall adjacent to a 
fallout field. 

It is instructive to discuss once niore the approx- 
imations of the preceding Sections, using the 
functions g,(6,¢) and g,(0). In figure 42.1, the 
points within the rectangle correspond to (0,ф) 
combinations for which g,(6,¢)=1; while the 
points outside the rectangle correspond to 0.(6,2) = 
0. The circles sunilarly divide the points in the 
plane into two sets, one giving g,(6,¢)=1, the 
other giving g;(@,¢)=0. But that of sketch (a) 


шы 
leads to a function g|(0) —5- | 9.(д,Ф)ф which 
0 


looks like Ше solid line of figure 42.1, (а), right; 
while the circular approximation (dashed line) 
resembles a step function. Note that the area 
under both curves must be about the same, since 
it corresponds approximately to the solid angle 
fraction. The circular sector and annular sector 
approximations are sunilarly compared with the 
functions they represent in the sketches (b) and 


(c), right, of the same figure. Thus all three 
approximations are seen to be combinations of 
step functions. 

Presumably greater accuracy 15 obtainable Бу 
use of approximations which involve more terms; 
but the greater effort is a serious impediment, and 
other sources of error tend to be large enough to 
make such refinenients have questionable sig- 
nificance. 
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Appendix A. Computations 


46. Preparation of Data: Use of Digital 
Computers 


To obtain most of the data presented іп this 
manuscript, 16 was necessary to rely on digital 
computer programs. А family of these programs 
was written for the IBM 704 computer at the 
National Bureau of Standards, using the SAP 
system of orders. Because so much of the desired 
data involved angular distributions of a fairly 
general type, many of these coniputer prograins 
involved extensive development and experimen- 
tation. Figure Al is a block diagram showing the 
system 1р its final form, and indicating by means 
of arrows the order followed. 

The first to be constructed was the general 
MOMENTS routine, whose output is а large 
tabulation of values for coefficients Z} (Е), as 
well as dose integrals 1005, for a given monoener- 
getic gamma ray source in a given material. Неге, 
Е is the 7th scattered energy and Ше indices л, 
[, ly refer to the n’th spatial moment, the РЪБ 
Legendre coefficient of the directional distribution, 
resulting from a plane source emitting gamma rays 
having initial directional distribution P, (cos 0), 
where 6, 1s the initial obliquity angle relative to 
the normal to the source plane [12]. The calcu- 
lation 1s carried out using photon wavelength (in 
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Compton units) as a variable. For most calcula- 
tions the basic integration interval was A\=0.025. 
This interval was doubled at АғА6--2, and re- 
doubled at Хм 4, to keep the total down to 
manageable proportions. For 1<[<12, values 
of n up to 12 were calculated; while for /,—0, Ше 
calculation was extended to n=18. Output was 


MOMENTS 
ABSTRACTOR 


ISOTROPIC PARAMETERS SLANT PARAMETERS 
ISOTROPIC DISTRIBUTIONS SLANT DISTRIBUTIONS 


Figure Al. Flow diagram for the computer programs. 


to magnetic tape, and totaled the order of 105 
numbers. After a number of time-saving modi- 
fications of the code, calculations for a single 
source energy and material took about one half 
hour of machine time. 

The purpose of the ABSTRACTOR code was 
to make a selection of data from the tape generated 
by the MOMENTS code. Two types of selection 
were considered, one essentially linear in the 
photon wavelength and the other essentially log- 
arithmic and including all the source wavelengths 
to be used. Both selections included the integral 
data, treating it as if it corresponded to spectral 
energies. (Four types of integral data were in- 
cluded: Data for the total dose, the total scattered 
dose, the once-scattered dose, and the multiply 
scattered dose.) By selecting about 1/10 of the 
scattered energies, it was possible to put on a 
single magnetic tape results for the following source 
wavelengths: 


м=0.05, 0.075, 0.1, 0.125, 0.15, 0.2, 0.25, 
0.3, 0.4, 0.5, 0.6, 0.8, 1.0, 1.2, 1.6, 2.0, 
24 3.0, 40 50 (0 90 112,00 150. 


Оп another tape was put data for the source wave- 
lengths 0.384, 0.434, and 0.77, corresponding to 
Соб and Сз!37. Initial calculations were for pene- 
tration in H;O. "These were followed eventually 
by calculations for concrete, and still later by 
calculations for Fe. 

Following the abstraction, the SORT code re- 
ordered the data from 


(E (lb n (КЕ))))) 
(Fo (lo (ELE; (n))))), 


because 1t. was thought that the tape search time 
would be cut down by having all the spatial 
moments In a group together. (This was probably 
an ineffective modification, all things considered.) 
From the basic moment-coefficient data it was 
possible to construct spatial and directional dis- 
tributions corresponding to many different source 
configurations. Of these configurations, two types 
had special importance: (a) those corresponding 
to isotropic sources (/,—0), and (b) those corre- 
sponding to plane sources einitting at а fixed 
obliquity angle 6, relative to the normal to the 
source plane. Accordingly a tape search sub- 
routine called SEEK was written, which, for 
specified №, Е, | would yield a set of moments 
describing one or the other type of distribution. 
Two codes were then constructed to determine 
actual distribution parameters. Inthe first, called 
ISOTROPIC PARAMETERS, the isotropic dis- 
tribution moments were used to calculate the 
coefficients of a power series representation of the 
point isotropic build-up factor. In other words, 
the (plane source) moments were turned into point 
source moments, and were used to determine 


to 
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coefficients A; in the representation 


] imax | 

ТЕ) n 24 А ;(ugr)!e Tinn (A1) 
of the Pth Legendre coefficient of the point 150- 
tropic source angular distribution. The А, co- 
efficients, for a specified set of „ЕЛ, апа Еу values, 
were punched on an IBM card, together with the 
descriptive 1,E;, апа № values. It may be of 
interest to note that this calculation involved a 
matrix inversion, which was accomplished in a 
very simple way by hand, using the bi-orthogonal 
polynomial system UF of Ref. 12. Equation 
(Al) refers also to the integral data, which, in 
fact, was utilized almost exclusively. 

In the second program of this type, called 
SLANT PARAMETERS, calculations were made 
of the constants 4°, 40, Bš, and В in a repre- 
sentation 


пе -Вег mar — В9г 

IP^(E,, cos 6, 2)=1/2 (35 де ПЗ Ae $ 
j=0 7=0 

(А2) 


of the lth harmonic coefficient of the angular 
distribution corresponding to a plane oblique (or 
slant) source at specified initial obliquity angle 00. 
In this calculation the values of Ве, B? were 
fixed at the value (ш/сов 06) and the values of 

° Во could be assigned arbitrarily. The other 
constants were all determined so that the moments 
of the approximate distribution agreed with the 
known values and with the known value of the 
source strength. The coefficients for “even” and 
"odd" component distributions, as indicated by 
the superscripts, were calculated separately and 
punched on separate IBM cards, together with 
identifying data for Ко, Е;, cos б, and a parity 
indicator which was zero or one. 

The calculation of slant distribution parameters 
was a fairly tricky business, because the B,’s were 
determined as the roots of polynomial equations. 
These roots could always turn out to be imaginary, 
signifying an oscillatory and therefore nonphysical 
distribution. The arbitrary B, coefficient was 
sometimes given a variety of values in the search 
for a reasonable, real set of roots. 

Many sets of parameter cards were punched, 
mostly corresponding to integral (dose) data, with 
1—0. Sets were prepared for values of l up to 
|--7 for isotropic sources, dose data, and various 
Е, spectral energies. Values for | greater than 
zero were not calculated for oblique sources, nor 
were calculations made for more than a very few 
Е, values. This was partly because most relevant 
shielding data did not require these details, and 
partly because the slant parameter calculations 
were expensive as well as tricky, and required care- 
ful examination of results. 

Tne final machine programs calculated spatial 
distributions, and in the case of the isotropic 
source cases, directional distributions. Specifi- 


cally, the SLANT DISTRIBUTIONS program 


could calculate penetration data for an arbitrary 
combination of source obliquities, source energies, 
and scattered energies, so long as these obliquities 
and energies were represented in the deck of 
binary parameter cards. In general, dose param- 
eter cards for the source energies already men- 
tioned were available, and the pow directions 
included cos b= + 1, "a Vo 2 2 . meli 

The ISOTROPIC DISTRIBUTIONS program 
was much more complicated. Initially, this pro- 
eram determined the total intensity or an angular 
distribution resulting from a Legendre sum, for a 
point isotropic source or a plane isoti opic source. 
Superposition of different source energies Was pos- 
sible. All the angular distributions for the fission 
source were obtained with this early form of the 
code. The resulting angular distributions were 
mostly not directly usable because of large oscilla- 
tions; it was a case of too few coefficients to 
describe a peaked distmbution. These angular 
distributions were smoothed by adding angular 
functions of the forms 


a(cosh B—cos دو‎ (I+-1/2)a,P,(cos 0), 


Point Isotropic Source, (A3) 
and 

7 
«(сов 0) 1g 818 us (14-1/2)6,Р(сов 0), 

=0 

Plane Isotropic Source. (А4) 


where а, В are arbitrary constants, and the a,’s 
and .b,’s are the expansion coefficients of the 
accompanying function. Itwill be seen that addi- 
tions of this type do not affect the values of the 
first 8, machine-calculated Legendre coefficients. 

The resulting, hand-calculated angular distri- 
butions are those of figures 26.1 and 26.3. Ве- 
cause of the work involved in obtaining these dis- 
tributions from the machine output, it appeared 
highly desirable to include methods for obtaining 
such distributions in the machine code itself. Ас- 
cordingly, a series of modifications was incor- 
porated in the code: (1) For the plane, isotropic 
calculations, the unscattered component was cal- 
culated exactly, and the scattered component was 
represented by a simple Legendre sum. (2) For 
the point isotropic calculations, the Legendre co- 
efficients were first calculated. Then the 1=3 and 
[—5 coefficients were used to obtain coefficients 
a and В for a function a(cosh 8— cos 0)??? so as to 
make the corresponding coefficients in the Legen- 
dre sum of that function equal to the known co- 
efficients. Finally, the angular distribution was 
calculated using the expression 


П (cos 0)—a (cosh B—cos 0) —3/2 


к; > (I+1/2)[D,—a)]P; (cos 8). (А5) 


(3) Fi ovisions were made to extend all саје ul: ations 
to include ring source, circular sourec, and circular 
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clearing source types, and another modification 
was made to provide for the determination of both 
unscattered and scattered components of the point 
isotropic source (total dose). 

Since this code provided for superposition of dif- 
ferent source energies, and for modification of the 
output angular distribution (Legendre sum only) 
to correspond to different reference axes, it made 
possible calculations for different spectra as well as 
the calculations of angular distributions incident 
on vertical walls (See eq 27.7). The angular dis- 
tributions in Appendix B were calculated using 
this code. 

Because of the extensive development work 
that accompanied the writing of these codes, they 
included many features which turned out not to 
be useful or desirable. They have never been put 
in such a form that general circulation appeared 
desirable. An effort is presently underway to 
rewrite the best features of this series, using the 
Fortran system, for general use. 


47. Preparation of Data: Hand Computations 


For the most part, the hand computations were 
perfectly straightforward integrations over the 
machine-produced data. One exception is to be 
noted: Integration of expressions eqs (27.11), 
(27.12), and (27.16) was accomplished in the most 
obvious way only once, namely for the fission 
data. Thereafter it was noted that one could 
obtain the same result in a far easier way by using 
integral data such as Ла, За and P, The ad- 
vantage 15 twofold: The azimuthal integration is 
carried out analytically and the obliquity inte- 
gration does not involve extremely peaked func- 
tions. The trick involves using an identity such 


as 
deos0l(d, cos0) = S(d)dS,,cos0 x 0, (A6) 


to write the expression for W, in the form 


=] 
/2о— =? |’ 
(А7) 


+ 





W,(d,1—cos а) == | dS, cos! 


Sa (d „1- sin a) 


where dS,=d|S,(d,w)]. 

In future calculations it would be highly dc- 
sirable to arrange it so that all the integrations 
әле carried out by machine. This will probably 
involve angular distribution output either in the 
form of punched cards or on magnetic tape. 


48. Comments on Accuracy 


In discussing the material of this manuscript in 
regards to accuracy, it is convenient to distinguish 
three types of error: (1) errors due to inaccurate 
input information; (2) errors duc to the procedure 
used to calculate distributions and produce the 
graphs; and (3) errors due to procedures for using 
the data, which are traceable to the inadequacies 
of the sehematization. 


It is clear that the use of a particular approxi- 
mate fission spectrum, namely that for 1.12 hour 
old fission products, to represent all practical situ- 
ations, is by far the largest error of the first type. 
The magnitude of such errors can be estimated 
roughly from data such as that in figure 26.6. 
For planning purposes the shielding factors may 
be defined relative to a particular spectrum. This 
doesn’t alter the fact that the data may not rep- 
resent an experimental situation exactly, but it 
does encourage a shift of attention to other parts 
of the problem, since then the validity of the 
spectrum and the accuracy of shielding factor cal- 
culations are considered separately. 


The other major type of input information 15 
the cross sectional data. We expect small inac- 
curacies in the differential scattering cross section 
due, e.g., to atomic binding corrections in the 
Klein-Nishina expression, to have negligible effect 
on the accuracy of the calculations. On the other 
hand, errors in the total cross sections which, gen- 
erally speaking, are expected to be less than 1 
percent in magnitude, can result in errors in the 
barrier factors at the deeper penetrations which 
are up to 10 times larger, i.e., 5 to 10 percent. 
A change in air density has the effect of a change 
In total cross section; and this can be brought 
about by temperature ог altitude variations. Fi- 
nally, not all data was prepared using precisely 
the conversion factor of eq 26.1; there were dif- 
ferences of a few tenths of a percent. For this 
reason, and due also to the use of overlays in 
preparing the graphs, the data are not expected 
to be internally consistent to better than about 
5 percent. 


Generally speaking, the attenuation data repre- 
sented m the barrier factors can be calculated 
with much greater accuracy than the angular dis- 
tributions which are used to obtain geometry 
factors. Calculations of barrier factors for point 
and plane isotropic sources, i.e., the L(X) and 
P(X) curves, are probably accurate to within 5 
percent, assuming that the input data is exactly 
correct. Considering the different sources of error 
together except for the spectrum assumption, these 
curves should be correct to within, say 8 percent. 
On the other hand, the S(d), S’(X), and W(.X,d) 
are appreciably less accurate; and it is difficult to 
assess quantitatively the accuracy of these curves 
even though some of the major sources of error 
are known. Perhaps for guidance purposes It may 
be useful to give "educated guesses” that excluding 
the assumption of spectrum, errors m S(d) niay 
be in the range of 10 to 15 percent; errors in 
W(X,d) for small X are perhaps 10 percent and 
for large .X may be two or three times larger; and 
errors in S (X) can be expected to increase roughly 
linearly with А, from 0 at S’=1 to perhaps as 
much as 100 percent at S’=107*.* 


9% Errors in S(d) stem from moderately rough calculations of the angular 
distribution; errors in W'(.V,d) come largely from the spectrum and anguler 
distribution assumed as source for the wall penetration; and errors in 85% А) 
are dua to its calculation as a small difference between large numbers, 


82 


Turning to angular distributions and geometry 
factors, the inherent inaccuracies of the calcula- 
tions can be seen visually in the Соб and Cs!” 
angular distributions in Appendix B, which have 
not been smoothed. The distributions of figures. 
26.1 and 26.3 originally looked very much the 
same but were arbitrarily smoothed with french 
curves. Errors in the range of 20 to 25 percent, 
perhaps even greater, can be expected at places 
in these data. However, the integrations which 
yield geometry factors improve the situation some- 
what. Generally speaking, we expect the geom- 
etry factors to be quite accurate when they are 
large and to approach the accuracy of the angular 
distributions from which they were obtained when 


they are very small. Thus, from right to left the . 


geometry factor curves vary from no error to 
errors of perhaps as much as 25 percent, with 
the L, and L, curves being more reliable than 
this because they were not obtained from angular 
distributions. Combinations of geometry and bar- 
rler factors, such as in figures 28.18, 28.19, and 
28.20 have different accuracy in different regions 
as determined by the separate errors inherent in 
the two ingredient factors. 

Finally, in regards to the utilization of the 
data, there are several comments which can be 
made. Use of water data to represent concrete 
gives results normally within the accuracy of the 
calculations if the 0.9 factor differences are cor- 
rectly applied. But it should be remembered that 
near the source, where the total contribution of 
scattered photons is a relatively small part of 
the detector response, the difference between water 
and concrete is negligible, certainly not as large 
as 10 percent ; and at great penetrations the differ- 
ence can be larger than 10 percent. Errors due 
to estimation of boundary effects probably ap- 
proach the magnitude of the total boundary effect. 
And errors due to the use of circular source data 
to represent rectangular source geometry factors 
wil vary greatly with the situation. They can 
perhaps be estimated by the use of several calcu- 
lations using different superpositions of data, as 
was indicated in the problem worked out In 
Section 43. 


49. Alternative Data Types and Procedures 


'The use of circular source data for computa- 
tion of geometry factors is only one of several 
possible procedures involving different data types. 
For example, in Section 44 it was indicated that 
geometry factors for rectangular source types 
could be calculated directly. 

An interesting alternative procedure has been 
explored by Berger, Hubbell, and Lamkin [29, 
30,41]. This consists of using the Legendre соећ- 
cients of the directional distributions rather than 
the directional distributions theniselves. If, m 
eqs (16.2) and (16.3) the dose angular distribu- 
tion. (D) is written in the form of a Legendre 
series, the detector response Ds takes the form 
of а sum of products m which one factor 1s a 





Legendre coefficient of the series and the other 
factor is the corresponding coefficient of the 
Legendre series representation of gs. Hubbell, 
Lamkin, and Bach [41] have prepared extensive 
tables of the gs Legendre coefficients to make 
such calculations feasible; but tables of Legendre 
coefficients for the relevant dose angular distribu- 
tions have never been prepared. 


The Legendre coefficient approach needs more 
attention. It is easy to apply, and has as its 
main limitation a tendency for the series of prod- 
ucts to converge slowly when the angular distri- 
bution has a sharp peak, as it often does (figs. 
26.1 and 26.3), and when the solid angle fraction 
of gs 1s at the same time small, as it often is. 
This limitation may prove moie apparent than real. 


Appendix B. Data for СО" and CS!” 


50. Introductory Comments 


Procedures for the calculation of the data rep- 
resented in the figures B1 through B50 have 
already been described. The main differences 
from the fission source data stem from the use 
of desk computers in that case for some calcula- 
tions performed by digital computer in the case 
of Соб and Сз!'37, the procedures involved not 
being entirely identical. 
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The figures given here are in the same sequence 
as the corresponding fission source figures in the 
main text; and data for the two source types have 
been placed adjacent to one another. 

No attempt has been made to analyze the obvi- 
ous similarities between the figures for the differ- 
ent source types. Jt is clear that many of the 
geometry factors are quite insensitive to source 
spectrum. 
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FIGURE Bl. Со-60, concrete. 
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FIGURE B2. Cs-187 
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IGURE B3. 
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FIGURE B7. Со-60, concrete. 
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Figure B8. Cs-187, concrete. 
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FIGURE B9. Со-60, backscatter from concrete. 
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Figure 812. Сз-187, H20. 
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FIGURE B17. Co-60, concrete. 
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FicunE B21. Со-60, concrete. 
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FIGURE B22. С5–187, concrete. 
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URE B23. Co-60, concrete. 
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FIGURE B25. Co-60, H2O data for l(d,cos0), concrete data for s(X,cos0,). 
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Ғтасве B26. Cs-187, Н2О data for l(d,cos0), concrete data for s( X,cos0,). 
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FIGURE B27. Со-60, concrete. 
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Figure B30. Cs-187, concrete. 
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FIGURE B32. Cs—137, concrete. 
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FIGURE B33. Co-60, H20. 
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FIGURE B36. Cs-137, H20. 
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FIGURE B39. Co-60, H20. 
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Figure B42. Cs-137, H20. 
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Figure B43. Со-60, concrete. 
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Figure B45. Со-60, concrete. 
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Figure B47. Со-60, concrete. 
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